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Space Telescope (HST) and have used them to trace the history of cosmic 
expansion over the last 10 billion years. These objects, which include 13 
spectroscopically confirmed SNe la at z > 1, were discovered during 14 epochs 



O ■ of reimaging of the GOODS fields North and South over two years with 

■ the Advanced Camera for Surveys on HST. Together with a recalibration of 

I ' our previous i/ST-discovered SNe la, the full sample of 23 SNe la at 2; > 1 

^ • provides the highest-redshift sample known. Combined with previous SN la 

^ ■ datasets, we measured H{z) at discrete, uncorrelated epochs, reducing the 
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uncertainty of H{z > 1) from 50% to under 20%, strengthening the evidence for 
a cosmic jerk-the transition from deceleration in the past to acceleration in the 
present. The unique leverage of the HST high-redshift SNe la provides the first 
meaningful constraint on the dark energy equation-of-state parameter at z > 1. 
The result remains consistent with a cosmological constant {w{z) = — 1), and 
rules out rapidly evolving dark energy {dw/dz >> 1). The defining property of 
dark energy, its negative pressure, appears to be present at z > 1, in the epoch 
preceding acceleration, with ~ 98% confidence in our primary fit. Moreover, the 
z > 1 sample-averaged spectral energy distribution is consistent with that of the 
typical SN la over the last 10 Gyr, indicating that any spectral evolution of the 
properties of SNe la with redshift is still below our detection threshold. 

subject headings: galaxies: distances and redshifts — cosmology: observations — 
cosmology: distance scale — supernovae: general 

1. Introduction 

The accelerating cosmic expansion first inferred from observations of distant type 
la supernovae (SNe la; Ricss ct al. 1998; Pcrlmutter et al. 1999) indicates unexpected 
gravitational physics, frequently attributed to the dominating presence of a "dark energy" 
with negative pressure. Increasingly incisive samples of SNe la at 2; < 1 have reinforced 
the significance of this result (Tonry et al. 2003; Knop et al. 2003; Barris et al. 2004; 
Conley et al. 2006; Astier et al. 2006). Using the new Advanced Camera for Surveys (ACS) 
and refurbished NICMOS camera on the Hubble Space Telescope (HST), our collaboration 
secured observations of a sample of the most-distant known SNe la. These half-dozen SNe 
la, all at 2; > 1.25, helped confirm the reality of cosmic acceleration by delineating the 
transition from preceding cosmic deceleration during the matter-dominated phase and by 
ruling out simple sources of astrophysical dimming (Riess et al. 2004b, hereafter R04). The 
expanded sample of 23 SNe la at 2; > 1 presented here are now used to begin characterizing 
the early behavior of dark energy. 

Other studies independent of SNe la now strongly favor something like dark energy 
as the dominant component in the mass-energy budget of the Universe. Perhaps most 
convincingly, observations of large-scale structure and the cosmic microwave background 
radiation provide indirect evidence for a dark-energy component (e.g., Spergel et al. 2006). 
Measurements of the integrated Sachs- Wolfe effect (e.g., Afshordi, Loh, & Strauss 2004; 
Boughn & Crittenden 2004; Fosalba et al. 2003; Nolta et al. 2004; Scranton et al. 2005) 
more directly suggest the presence of dark energy with a negative pressure. Additional, 
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albeit more tentative, evidence is provided by observations of X-ray clusters (Allen et al. 
2004) and baryon oscillations (e.g., Eisenstein et al. 2005). 

The unexplained existence of a dominant, dark-energy-like phenomenon presents a 
stiff challenge to the standard model of cosmology and particle physics. The apparent 
acceleration may result from exotic physics such as the repulsive gravity predicted for a 
medium with negative pressure or from entirely new physics. The explanation of strongest 
pedigree is Einstein's famous "cosmological constant" A (i.e., vacuum energy; Einstein 
1917), followed by a decaying scalar field similar to that already invoked for many infiation 
models (i.e., quintessence -Wcttcrich 1995, Caldwell, Dave, & Steinhardt 1998; Peebles 
& Ratra 2003). Competitors include the Chaplygin gas (Bento, Bertolami, & Sen 2002), 
topological defects, and a massless scalar field at low temperature. Alternatively, alterations 
to General Relativity may be required as occurs from the higher-dimensional transport of 
gravitons in string theory models (Deffayet et al. 2002) and braneworlds, or by finely-tuned, 
long-range modifications (e.g., Cardassian type, Frccsc 2005; or Carroll et al. 2004; see 
Szydlowski, Kurek, & Krawiec 2006 for a review). Empirical clues are critical for testing 
hypotheses and narrowing the allowed range of possible models. 

SNe la remain one of our best tools for unraveling the properties of dark energy because 
their individual measurement precision is unparalleled and they are readily attainable in 
sample sizes of order 10^, statistically sufficient to measure dark-energy-induced changes to 
the expansion rate of ~1%. Specifically, the equation-of-state parameter of dark energy, w 
(where P = wpc^) determines both the evolution of the density of dark energy. 



where Pde,o is the present dark-energy density. Measuring changes in the scale factor, a, 
with time from the distance and redshift measurements of SNe la, 

di{z) [to df dz' 



constrains the behavior of w{a) or w{z) and is most easily accomplished at 2; < 2 during the 
epoch of dark-energy dominance. 

Ideally, we seek to extract the function w{z) for dark energy or its mean value at a 
wide range of epochs. Alternatively, we might constrain its recent value Wq = w{z — 0) 
and a derivative, dw / dz = w' , which are exactly specified for a cosmological constant to be 
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and its gravitational effect on expansion. 
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(—1,0). Most other models make less precise predictions. For example, the presence of a 
"tracker" dark-energy field whose evolution is coupled to the (decreasing) dark matter or 
radiation density may be detected by a measured value of > or ^(2; > 1) ~ 0. In truth, 
we know almost nothing of what to expect for w{z), so the safest approach is to assume 
nothing and measure w{z) across the redshift range of interest. SN la at 2; > 1 are crucial 
to constrain variations of w with redshift. These measurements can only be made from 
space, and we report here on that endeavor. We have discovered and measured 21 new SN 
la with HST and used them to constrain the properties of the dark energy. We present the 
follow-up spectroscopy and photometry of the new SNe la in §2, light-curve analysis and 
cosmological constraints in §3, a discussion in §4, and a summary in §5. 



2. Further Discoveries and Data Reprocessing 

2.1. ACS Searches, 2003-2005 

In HST Cycle 11 (July 2002 - June 2003) we initiated the first space-based program 
designed to find and monitor SNe (R04). Our search was conducted by imaging the two 
high-latitude fields of the GOODS Treasury Program (the Chandra Deep Field South and 
the Hubble Deep Field North) with 15 ACS pointings, 5 times each at 6 or 7-week intervals 
(chosen to match the rise time of SNe la at 2; ~ 1). Multiple exposures in the F850LP 
bandpass were differenced to find transients. Contemporaneous color measurements and 
host-galaxy photometric redshifts were utilized to identify promising candidate SNe la at 
2; > 1 for target-of-opportunity (ToO) follow-up observations (Riess et al. 2004a). 

In HST Cycles 12 and 13 (July 2003 - June 2005) we continued our past efforts, 
imaging the GOODS fields (Giavalisco et al. 2004, Strolger et al. 2004) 14 more times and 
following newly discovered SN la candidates with HST. The observational methods we used 
in Cycles 12 and 13 were very similar to those of Cycle 11 and are extensively described by 
R04. Readers are directed to that publication for the sake of brevity; here we describe only 
changes to our observing and candidate-selection strategies. 

To improve our search efficiency with finite observing time, we compressed our primary 
search exposures from 2 orbits to 1. We omitted the short- wavelength filters, whose use in 
the GOODS program was primarily to constrain the properties of galaxy formation. Our 
typical search sequence contained four dithered exposures of length 400 s in F850LP and 
one 400 s exposure in F775W for a typical orbit instead of the previous four exposures of 
500 s length in F850LP (as well as two such exposures in each F606W and F775W). At 
times when we had unusually long orbits, especially approaching the continuous viewing 
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zone for the HDFN, we included a ~200 s exposure in F6061^ as well to help constrain SN 
type and redshift. 

Although our net exposure time was reduced by 20% from the Cycle 11 observations, 
the accumulation of additional template images for the host galaxies without the supernova 
during the GOODS program increased their total depth and more than compensated for 
the slightly decreased signal in this search. 

We benefited from two modest improvements in our ability to identify high-redshift 
SNe la before initiating a ToO over the strategy utihzed by R04. The first was the 
availability of spectroscopic host-galaxy redshifts obtained before the appearance of a SN 
candidate (e.g., from our Keck and other spectra assembled for the GOODS catalog, and 
from "Team Keck"; Wirth et al. 2004). The second was the result of building up a 3 year 
baseline of variability for each galaxy. This allowed us to distinguish between transient 
signals caused by supernovae and transient noise introduced by the variable emission from 
active galactic nuclei. 

The 21 new SNe la are listed in Table 1 and their discovery images are shown in Figure 
1. Internal names for the SNe are used in the following. Color images of all of the SN la 
host galaxies (where filters F850LP, F606W, and FiSbW correspond to red, green, and 
blue, respectively) are shown in Figure 2, with the position of the SN indicated in each case. 

2.2. Photometry 

Our follow-up observations of candidate SNe la were similar to those previously 
obtained and described by R04. After the search phase, all images were reprocessed using 
up-to-date reference files and the CALACS pipeline in the STSDAS package in IRAF.^ 
Because improvements have been made to the standard CALACS pipeline (including 
an improved distortion map and more precise flatfield images), we also reprocessed all 
preceding SN data from our first year. 

To produce light curves, we developed an automated pipeline to retrieve all images 
obtained by ACS and NICMOS in the GOODS fields from the Multimission Archive 
(reprocessing all data from R04 as well as all of the new frames). The enhanced set of 
images available from before and after each SN increased the depth of the constructed 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Table 1. SN Discovery Data 



Name 


Nickname^ 


JD-2.400,000 


SN a(J2000) 


SN (5(J2000) 





HST04Sas 




53148.7 


12:36:54.11 


+62-08-22 76 


HST04Man 


Manipogo 


53146.3 


12:36:34.81 


+62:15:49.06 


HST04YOW 


Yowie 


53145.2 


12:36:34.33 


+62:12:12.95 


HST04Pat 


Patiixent 


53196.3 


12:38:09.00 


+62:18:47.24 


HST04Tha 


Thames 


53196.5 


12:36:55.17 


+62:13:04.05 


HST04Cay 


Cayanne* 


53245.2 


12:37:27.11 


+62:12:07.68 


HST04Eag 


Eagle 


53284.9 


12:37:20.75 


+62:13:41.50 


HST04Haw 


Hawk 


53284.0 


12:35:41.16 


+62:11:37.19 


HSTOSSpo 


Spock 


53376.8 


12:37:06.53 


+62:15:11.70 


HST04Mcg 


McGuire 


53265.0 


03:32:10.02 


-27:49:49.98 


HST04Gre 


Greenberg 


53265.0 


03:32:21.49 


-27:46:58.30 


HST04Omb 


Ombo 


53356.3 


03:32:25.34 


-27:45:03.01 


HST04Rak 


Rakke 


53356.3 


03:32:18.15 


-27:44:10.55 


HST04Kur 


Kurage 


53355.2 


03:32:36.03 


-27:51:17.66 


HST05Dic 


Dickinson 


53472.5 


12:35:49.61 


+62:10:11.96 


HST05Fer 


Ferguson 


53472.6 


12:36:25.10 


+62:15:23.84 


HST05Koe 


Koekemoer 


53472.7 


12:36:22.92 


+62:18:23.20 


HSTOSStr 


Strolger 


53474.5 


12:36:20.63 


+62:10:50.58 


HST05Gab 


Gabi 


53474.5 


12:36:13.83 


+62:12:07.56 


HST05Red 


Rcdford 


53426.6 


12:37:01.70 


+62:12:23.98 


HST05Lan 


Lancaster 


53427.1 


12:36:56.72 


+62:12:53.33 


HS05Zwi 


Zwicky 


53403.2 


03:32:45.65 


-27:44:24.30 



*not positively identified as SN la 
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template image (without SN light) for each SN, thereby increasing the subsequent precision 
of the SN photometry. 

Two modest improvements to measurements of SN flux were made from the methods 
described by R04. First, due to backside scattering from the ACS WFC CCD mounting, the 
point-spread function (PSF) redward of ~ 8500 A displays a halo of scattered hght which 
needs to be quantified for precise photometry in the F850LP bandpass. R04 used a red star 
with i — z color similar to a SN la at peak at z ^ 1.3 {i — z = 1.0 mag, Vega) as a model PSF 
for SN photometry. However, variations in SN rcdshift and phase around that model make 
the match with a single PSF at a single color inexact resulting in photometric errors ranging 
from ±0.05 mag for our sample of phases and redshifts. Here, we adopted an individualized 
method for the PSF modeling of high-redshift SNe la. As described in detail by Casertano 
& Riess (2007), we used 6 observed monochromatic PSFs from ACS narrow-band filters 
between 8150 A and 10500 A to derive an empirical wavelength dependent correction to 
the average PSF model. At these monochromatic endpoints the difference between (nearly) 
infinite-aperture photometry and that determined from a fit to the previously used red-star 
PSF is —0.14 mag and +0.46 mag, respectively, with no difference occurring near 9300 A. 
These differences were evaluated for all monochromatic wavelengths. At every SN redshift 
and phase, the same representative model SN la spectral energy distribution (SED) used to 
derive K-corrections was used to integrate the monochromatic PSFs within the F850LP 
band response and determined their photometric differences. In Figure 3 we show the 
required photometric correction as a function of SN la age and redshift. Spatial variations 
in the ACS WFC PSF are small (Krist 2003) and variations in the position of the SNe 
result in a PSF noise of 1-2%, well below the sky and read noise. 

Second, we adopted an updated zeropoint for the HST ACS "Vegamag" photometric 
system from Sirianni et al. (2005) which uses the spectrophotometry of Vega from Bohlin 
& GiUiland (2004) and improved quantum efficiency data for ACS to set Vega to 0.00 
mag in all passbands. We then iitilizcd the same magnitudes of Vega assumed by Landolt 
(1992) to calibrate the HST photometry on the same Vega system as the nearby SN la 
sample. Thus our data is calibrated on the Landolt system. The resulting zeropoints 
corresponding to 1 electron s"^ are F850LP = 24.35, F775W = 25.28, and F606W = 26.43 
mag. These are fainter than those used in R04 by 0.02 mag. For NICMOS we use the values 
of FUOW = 22.92 mag and FIGOI^ = 22.11 mag given by the STScI NICMOS Handbook. 

An improvement in our ability to estimate photometric errors was provided by the 
increased set of host images long before or after the appearance of the SNe. For each such 
image we added and measured an ideal PSF at a range of magnitudes at the same pixel 
position as the site of the SN and derived statistics from the recovered magnitudes such as 
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their dispersion and their bias at a given magnitude.^^ 

The cahbration of HST NICMOS NIC2 photometry has changed since R04 due to 
the recent detection of an apparent non-hnearity in all NICMOS detectors. This effect 
was initially discovered by Bohlin, Lindler, & Riess (2005) and has now been cahbrated to 
good precision by de Jong et al. (2006) using pairs of star-cluster images obtained with the 
flatfield lamps switched on and off. The result for NIC2 is a reduction in the apparent flux 
by 0.06 and 0.03 mag per dex in FllQW and Fl&QW , respectively.^^ 

For faint sources such as high-redshift SNe la, the nonlinearity in apparent flux pertains 
to the flux difference between the calibration stars (G191B2B and P330E, both ~12 mag in 

J and H) and the faint SN. For SN fluxes fainter than the sky level (all of those presented 
here) the nonlinearity pertains to the flux difference between the calibration stars and the 
sky level, below which any additional nonlinearity of sources is effectively "quenched". 
For SN la plus host ffuxes near or below the sky (a typical sky level is 0.17 electron s~^ 
in FllQW and 0.14 electron s^^ in FIQQW), the correction we calculate and apply is 
0.220 mag brighter (than the uncorrected zeropoints) in FllOl^ and 0.086 mag brighter in 
F1601^. Interestingly, the change in distance modulus from R04 due to these corrections is 
mitigated by their compensating effect in distance and reddening. 

SNe measured in our detection and follow-up observations are generally of high enough 
signal-to-noise ratio (S/N > 10) to support the use of magnitudes (where a magnitude is 
zeropoint — 2.5 logigflux) without skewing the interpretation of photometric errors. The 
infrequent observation (less than 5% of observations) with (S/N < 5) on the post- maximum 
light curve of any SN has insignificant weight when used in conjunction with the higher 



We found biases to exist only for the case of bright, sharp hosts which were "softened" by image 
interpolation in the registration of the template. The worst case hosts such as 2003XX and 2002hp seen in 
Figures 1 and 2 without correction would bias the magnitudes too bright by as much as 0.2 mag at the end 
of the observed light curve. We identified two valid solutions to this bias which gave consistent results: (1) 
equally soften the SN image by reinterpolation, or (2) estimate the bias from the simulated SNe la. This 
latter option is preferable as it is nearly noise-free. 

^^The correction is described by dc Jong et al. (2006) in the form: count rate a flux", where 
a = 1.025 ± 0.002 and a = 1.012 ± 0.006 for F110T4^ and F160VK, respectively. 

^■^Because the change in the bluer band [FllQW) is larger, the net change in the distance moduli from 
R04 to first order is approximately AFllOW^ - R{AFllQW - AF160W^), which for a reddening ratio i? » 3 
is 0.17 mag. In practice, a second-order change occurs as the apparent color impacts other parameters in 
the fit as well as the individual K-corrections. For the highest-redshift SNe which rely heavily on NICMOS 
data, the average change from R04 was found to be 5% (0.10 mag) closer in distance, but depends on the 
individual SN la. 



-9- 



signal-to-noise ratio points. However, our pre-discovery images require the use of flux 
measurement and error estimates due to their greater leverage on the determination of the 
time of maximum. For these images we measured the sky-subtracted flux (electrons per 
second) in a 0.17" radius aperture within which the zeropoint for a 1 electron per second 
source is is 24.68 mag for F850LP. 

Our flnal photometry for all SNe la is listed in Table 2 and shown compared to the 
individual multicolor light-curve shape (MLCS) flts in Figure 4. 

2.3. Spectroscopy and SN Identification 

The ACS grism spectroscopy we obtained for the SNe la is listed in Table 3 (they 
will also be made available at the University of Oklahoma supernova spectra database, 
SUSPECT). In general, our reduction and analysis methods were the same as in R04 with 
the following exceptions. A more realistic skyflat was used to separate the contributions 
to each pixel from the sky and the source as described in Pirzkal et al. (2005). We also 
utilized a Lanczos kernel in the drizzling procedure to decrease the effective size of the pixel 
convolved with the images and improve the separation for cases where the SN and host 
were very close together (e.g., SN HST04Sas). 

The spectra we obtained and used to classify the SNe are shown in Figure 5. As in 
R04, to classify the SNe, the detected SN spectra were cross-correlated with template 
spectra (after removal of the continuum) to identify their type and redshift using the 
"SNID" algorithm (Tonry et al. 2003, Blondin et al. 2006). For the cases listed in Table 3 
for which narrow-line host emission was identifled, the redshift was constrained to the value 
determined from the host emission before cross-correlation, improving the signiflcance of 
the cross-correlation peak. For all 12 spectra shown in Figure 5, SNID provided a signiflcant 
classification for each as type la. Although the diagnostic used by the SNID algorithm 
relies on the whole spectrum, the majority of these SNe can also be classified as type la 
from the presence of Si II absorption at 4130 A (Coil et al. 2000). Specifically, evidence of 
Si II absorption was seen in the two highest-redshift spectra presented here. Broad Ca II 
absorption near 3750 A is visible in all the spectra as well, but this feature is less secure 
than Si II for SN la classification due to its appearance in the spectra of SNe Ic (Fihppenko 
1997). For the highest redshift spectrum shown, HST04Sas, FllOI^, F160W, and F205W 
band NIC2 imaging of the host was obtained to constrain the phot-z of the red, elliptical 
host to z — 1.4 ± 0.15 which compares well with the SN value of 2; = 1.39 ± 0.01. 

For SN HST04Cay cross-correlation peaks exist, but none with high significance and 
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the redshift of the host is not known. For HSTOSRed, the spectral match to a SN la at the 
host ^ = 1.189 in Table 3 is fair, but not secure. SN 2003XX, identified by R04 as an SN la 
on the basis of its elliptical host, is now spectroscopically classified as an SN la due to the 
acquisition of an ACS grism spectroscopic galaxy-only template on 2005-03-07 to subtract 
from the grism spectrum of the SN obtained on 2003-04-16 and heavily contaminated by 
galaxy light. We also analyzed the ACS grism spectra of two more apparent high-redshift 
SNe observed by another collaboration (GO-9729; P.I. Perlmutter) during the course of 
Cycle 12 in the GOODS fields (SN150G a = 12:37:09.456, 5 = +62:22:15.59; SN150I at a = 
12:37:51.533, 5 — -1-62:17:08.24) but failed to find a significant peak in the cross-correlation 
so we cannot determine their type or redshift. 

Because of the large range in the quality and breadth of the photometric record of 
individual SNe la, R04 developed a two-tiered approach to the confidence of our SN la 
identifications. By distinguishing or selecting data based on objective criteria of their good 
quality we can mitigate systematic errors caused by undersampling of light curves and 
misidentifications without introducing cosmological biases. 

To summarize this approach: we classify as "high-confidence" (hereafter as "Gold") 
SNe la those with a compelling classification and whose photometric record is sufficient 
to yield a robust distance estimate easily characterized by its measurement uncertainty. 
"Likely but uncertain SNe la" (hereafter "Silver") are those with an aspect of the spectral 
or photometric record which is absent or suspect and whose distance error is described with 
a caveat rather than a quantitative uncertainty. As in R04, the three primary reasons for 
rejecting a SN la from the Gold set are that (1) the classification, though plausible, was 
not compelling, (2) the photometric record is too sparse to yield a robust distance (i.e., the 
number of model parameters is approximately equal to the effective number of samplings 
of the fight and color curves, and (3) the extinction is so large as to be uncertain due to 
our ignorance of extragalactic extinction laws. Although R04 set this extinction threshold 
to 1 mag of visual extinction, here we adopt the more conservative threshold of 0.5 mag 
utilized elsewhere (Miknaitis ct al. 2007; Tonry et al. 2003; Riess et al. 2005). SNe with 
two liens against its confidence (i.e., those rejected from the "Silver" set) are not included 
in the remainder of this paper (e.g., HST05Cay). 

The measured SN la distances for the full Gold and Silver set of i/ST-discovered SNe 
are given in Table 6, including the revised distance measurements of the HST objects first 
presented by R04. These distances use the same distance scale as in R04. To this sample 
we add the same ground-discovered sample of SNe la employed in R04 for the following 



- 11 - 



analyses ■'^ 

An important and recent addition to the ground-discovered sample is the first-year 
SNLS dataset from Astier et al. (2006) containing 73 new SNe la. We sought to define 
a similarly high-confidence subsample from this set using the information available from 
Astier et al. (2006) and Howell et al. (2005). 

We removed from the full-sample 2 events discarded by the SNLS (SNLS-03D4au and 
SNLS-03D4bc), as well as the 15 objects classified by Howell et al. (2005) and Astier et 
al. (2006) as "SN la*," meaning they are "probable SN la" but among the "least secure 
identification," and another type (such as a SN Ic) is not excluded. As discussed by Astier 
et al. (2006), most of the objects at ^ > 0.8 suffer from imprecise color measurements 
which dominate the distance error and would lead to a distance bias after application of 
a low-extinction cut. However, the low-extinction cut is valuable to reduce the sensitivity 
to unexpected or evolving extinction laws. Our solution is to first eliminate objects with 
highly uncertain color measurements, (Jcoior > 0.15 mag, which removes the following 6 
objects: SNLS-03Dlew, SNLS-03D4cn, SNLS-03D4cy, SNLS-04D3cp, SNLS-04D3dd, and 
SNLS-04D3ny. The distance uncertainty of these objects is typically 2 to 3 times that of the 
rest of the objects and thus the value of these 6 is very low (roughly the equivalent of losing 
a single well-measured SN la). Lastly, the same color cut of Ay < 0.5 mag used for all Gold 
sample SNe removes 3 objects: SNLS-03Dlgt, SNLS-03D3ba, and SNLS-04D2gc. Thus, 47 
high-confidence SNLS objects remain classified as "Gold". The SNLS SNe were fit with our 
present MLCS2k2 algorithm to estimate their distances. Alternatively, we found that after 
the addition of 0.19 mag to distances measured to the SNLS SNe by Astier et al. (2006) (to 
account for the arbitrary choice of distance scale and determined from low-redshift SNe in 
common), the agreement between our MLCS2k2 measured distances on the distance scale 
used in R04 and here and those from Astier et al. (2006) were consistent in the mean to 
better than 0.01 mag. Thus, either set of distances measured to the SNLS SNe provide a 
comparable and suitable addition to our cosmological sample. We verified (next section) 
that either yields the same inferences for w{z) and provide results fitting the SNLS SNe 
with either method in Table . 



^^As in R04, past SN data such as light curves and spectra from Perlmutter et al. (1999) and spectra 
from the sample of Knop et al. (2003) remain unavailable. Thus we resort to the same reliance as in R04 on 
their published distances normalized to a consistent distance scale using SNe in common and classification 
confidences from these authors. In §4 we also consider the impact of rejecting these and other older data. 

^^We made use of the Astier et al. (2006) distances for our primary fit as their light curves were not 
initially published by Astier et al. and only recently made available during the preparation of this work. 
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We have not made use of any of the seven high-redshift SNe la from Krisciunas et al. 
(2005) due to the apparently biased selection of the sample, as discovered by Krisciunas et 
al. but not yet fully modeled and corrected. 

Thus, a simple description of the full sample used here for the cosmological analyses 
consists of the addition of Table 5 from R04 to Table 6 provided here with the revised 

distances to the SNe in common (i.e., the leading 20 SNc in Table 6) superceding those 
given in R04 Table 5, and the SNLS objects from Astier et al. (2006) identified in §2.3. 
Either the Astier et al. (2006) provided distances (with the aforementioned addition of 
0.19 mag) or our own fits are comparable and suitable. The sample can also be found at 
http://braeburn.pha.jhu.edu/~ariess/R06 or upon request to ariess@stsci.edu. 

Upcoming revisions to the ground-discovered samples and improvements to the 
distance-fitting algorithms are expected and will change the membership and distance 
measures in the full cosmological sample and these should be considered before construction 
of a cosmological sample of SNe la. In response to such improvements we will attempt 
to provide distance estimates to the HST-discovered sample with updated fitting tools or 
distance scaling as warranted at http://braeburn.pha.jhu.edu/~ariess/R06 or upon request 
to ariess@stsci.edu. 



Table 2. SN la Imaging 



Date" Vega Mag Epoch(rest) K-Corr 



HST04Pat {z = 0.97) 





F775W 




F775W — *• U 


^ fj -L rj \J , 


23 86 fO 041 


-1.1 


-0 63^0 10) 


oo^yjo. -L 


.^O.OUIU.-LL/ i 


3 3 


U. 1 Ol U.UO ( 




FR^OLP 




_F850LP B 


53196 3 


23 77f0 04) 


-1.1 


-1 25f0 05) 




9^^ SI (Ci 04^ 


3 3 




53205 4 


23 77(0 04) 


3.4 


-1 28f0 04) 


53225 3 


24 45 fO 04) 


13.5 


-1 35f0 05) 


53377.7 


27.01(0.40) 


90.9 


-1.49(0.01) 




HST04Mce: (z 


= 1.37) 






J. 00\J-L/-L 




PS5nr p ^ TT 


53265 


24 M(n 04) 


5.9 


-1 01 (0 03) 




94 fil (Tl nfil 




-l.UU^U.U^ J 


53285 6 


25 39^0 14) 


14.6 


-1 1 1 fO 02) 


53294 7 


25 69f0 17) 


18.4 


-1 12f0 01) 


53306 5 


25 95('0 14) 


23.4 


-1 15f0 02) 

-L ft -L ^ 1 W ft V/ ^ f 


53312 3 


27 OOfO 50) 


25.9 


-1 16(0 02) 




i^llOVF 




FllOVF B 


53277 6 


24 22('0 05) 


11.2 


-1 71 ro 07) 






1 <^ n 

±0.\J 


- ± ft ( ( I u. -'-■^ i 




F160W 




F160VF R 


53277.7 


23 99f0 05) 


11.3 


-1 93f0 05) 


53285 6 


24 37f0 12) 


14.6 


-1 90f0 02) 


53217.9 


flux:0. 180(0.06) 


-13.9 


( ) 




xlo ± uorer \^Z 


— i.uz j 






F775W 




i^775W^ ^ t/ 


53486.8 


24.55(0.10) 


10.1 


-0.82(0.04) 




F850LP 




F850LP ^ B 


53472.6 


23.58(0.03) 


3.1 


-1.39(0.02) 


53486.8 


24.14(0.04) 


10.1 


-1.40(0.01) 


53493.8 


24.28(0.06) 


13.6 


-1.40(0.02) 


53504.3 


24.99(0.08) 


18.8 


-1.43(0.03) 


53513.2 


25.24(0.08) 


23.2 


-1.45(0.01) 


53527.3 


26.08(0.17) 


30.2 


-1.45(0.01) 




FllOW 




Fiiow^y 


53486.4 


23.73(0.05) 


9.9 


-1.40(0.05) 


53491.2 


23.87(0.06) 


12.3 


-1.37(0.05) 


53427.7 


flux:-0.03(0.03) 


-19.0 


(-) 




HST05Koe {z 


= 1.23) 






F775W 




J^775W" ^ U 
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Table 2 — Continued 



Date" 


Vega Mag 




Epoch(rest) 


K-Corr 


53472.7 


25.18(0.05) 




8.5 


-0.15(0.05) 


53486.8 


26.11(0.10) 




14.8 


-0.07(0.03) 




F850LP 






F850LP B 


53472.7 


24.34(0.06) 




8.5 


-1.43(0.01) 


53486.8 


25.11(0.07) 




14.8 


-1.40(0.04) 


53493.8 


25.34(0.11) 




18.0 


-1.34(0.08) 


53504.3 


25.90(0.15) 




22.7 


-1.24(0.04) 


53513.2 


26.10(0.16) 




26.7 


-1.19(0.02) 


53527.3 


27.10(0.40) 




33.0 


-1.18(0.03) 




mow 






FllQW 


53485.3 


24.48(0.07) 




14.1 


-1.46(0.08) 


53491.3 


24.66(0.09) 




16.8 


-1.42(0.09) 






HST05Dic {z -- 


= 0.638) 






F606W 






FQmw U 


53472.5 


23.42(0.02) 




4.3 


0.160(0.03) 




F775W 






F77bW -> B 


53472.5 


22.61(0.02) 




4.3 


-0.96(0.03) 




F850i:P 






F850LP V 


53472.5 


22.44(0.02) 




4.3 


-1.06(0.02) 


53530.8 


23.97(0.02) 




39.9 


-1.09(0.01) 






HST04Gre {z 


= 1.14) 






F775W 






F775W U 


53265.1 


23.61(0.04) 




1.0 


-0.46(0.01) 


53288.5 


24.58(0.05) 




12.0 


-0.46(0.01) 


53299.5 


25.22(0.07) 




17.1 


-0.46(0.01) 


53313.2 


25.57(0.10) 




23.5 


-0.47(0.01) 




F850LP 






P850LP ^ B 


53265.1 


23.25(0.02) 




1.0 


-1.51(0.02) 


53275.6 


23.38(0.02) 




6.0 


-1.49(0.02) 


53288.5 


24.00(0.05) 




12.0 


-1.47(0.03) 


53299.5 


24.68(0.10) 




17.1 


-1.43(0.05) 


53313.2 


24.98(0.07) 




23.5 


-1.36(0.02) 


53356.3 


26.26(0.31) 




43.7 


-1.32(0.03) 


53378.4 


26.31(0.38) 




54.0 


-1.37(0.02) 


53388.6 


25.94(0.26) 




58.8 


-1.36(0.01) 


53398.1 


26.79(0.32) 




63.2 


-1.37(0.01) 


53404.2 


26.10(0.23) 




66.0 


-1.37(0.01) 




FUOW 






FllOW V 


53290.5 


23.62(0.05) 




12.9 


-1.49(0.06) 


53301.2 


23.90(0.05) 




17.9 


-1.42(0.05) 


53217.2 


flux:-0.05(0.04) 


-21.2 


(-) 
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Table 2 — Continued 



Date" Vega Mag Epoch(rest) K-Corr 



HST04Omb {z = 0.975) 





F77f>W 






^^775^" ^ U 


53356.3 


23.82(0.05) 




-5.5 


-0.58(0.05) 


53373.7 


23.81 0.05) 




3.2 


-0.74(0.07) 


53378.5 


24.06(0.06) 




5.7 


-0.79(0.05) 


53388.6 


24.46(0.06) 




10.8 


-0.86(0.05) 


53398.1 


24.77(0.07) 




15.6 


-0.93(0.02) 


53404.1 


25.09(0.08) 




18.6 


-0.94(0.01) 










i^850LP B 


53356.3 


23.78(0.04) 




-5.5 


-1.26(0.02) 


53373.7 


23.58(0.03) 




3.2 


-1.30(0.03) 


53378.4 


23.61(0.03) 




5.6 


-1.32(0.02) 


53388.6 


23.99(0.05) 




10.8 


-1.34(0.03) 


53398.1 


24.31(0.06) 




15.6 


-1.39(0.05) 


53404.2 


24.40(0.06) 




18.7 


-1.44(0.06) 


53407.9 


24.50(0.07) 




20.6 


-1.47(0.05) 




Fiim 






FllQW 


53368.9 


23.31(0.04) 




0.8 


-1.43(0.02) 


53377.8 


23.54(0.04) 




5.3 


-1.41(0.06) 






HST05Red iz = 


1.19) 












F775W U 


53426.6 


25.52(0.05) 




12.0 


-0.21(0.03) 


53436.6 


25.92(0.07) 




16.6 


-0.19(0.01) 




F850LP 






F850LP B 


53426.6 


24.42(0.04) 




12.0 


-1.38(0.04) 


53436.6 


24.74(0.06) 




16.6 


-1.34(0.07) 


53439.6 


25.10(0.11) 




18.0 


-1.30(0.07) 


53445.5 


25.47(0.14) 




20.7 


-1.26(0.07) 


53454.8 


26.22(0.21) 




24.9 


-1.20(0.02) 


53463.3 


26.58(0.25) 




28.8 


-1.19(0.03) 


53471.7 


26.75(0.34) 




32.6 


-1.21(0.02) 


53472.9 


26.65(0.20) 




33.2 


-1.21(0.03) 




FllQW 






FllQW^V 


53440.4 


24.31(0.05) 




18.3 


-1.30(0.05) 




Fimw 






FIQQW R 


53440.5 


23.51(0.08) 




18.4 


-1.69(0.13) 






HST05Lan {z = 


1.23) 






i^775W^ 






FllbW U 


53427.1 


25.50(0.05) 




-2.4 


-0.24(0.05) 


53436.6 


25.75(0.05) 




1.8 


-0.19(0.05) 


53436.7 


25.88(0.10) 




1.8 


-0.19(0.05) 
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Table 2 — Continued 



Date" 


Vega Mag 




Epoch(rest) 


K-Corr 




F850LP 






F850LP B 


53427.2 


24.75(0.10) 




-2.4 


-1.40(0.03) 


53436.6 


24.56(0.08) 




1.8 


-1.37(0.03) 


53439.6 


24.72(0.09) 




3.1 


-1.36(0.04) 


53445.5 


24.96(0.06) 




5.8 


-1.34(0.02) 


53454.8 


25.29(0.08) 




9.9 


-1.33(0.02) 


53463.3 


25.69(0.18) 




13.7 


-1.31(0.04) 


53471.7 


25.96(0.26) 




17.5 


-1.26(0.08) 


53473.6 


26.78(0.42) 




18.4 


-1.24(0.08) 




FllOW 






FllOW 


53439.4 


24.40(0.06) 




3.0 


-1.56(0.05) 


53445.4 


24.43(0.06) 




5.7 


-1.53(0.06) 




F160W 






F160W R 


53438.4 


24.36(0.14) 




2.6 


-1.80(0.07) 


53444.4 


24.43(0.14) 




5.3 


-1.73(0.06) 






HST04Tha {z ■ 


= 0.954) 






F775W 






F775W U 


53196.5 


24.19(0.05) 




9.8 


-0.91(0.05) 


53207.7 


25.16(0.10) 




15.5 


-0.99(0.02) 




i^850LP 






F850LP B 


53196.5 


23.87(0.03) 




9.8 


-1.37(0.03) 


53207.7 


24.37(0.06) 




15.5 


-1.43(0.06) 


53216.3 


24.72(0.11) 




19.9 


-1.51(0.06) 


53221.7 


25.18(0.10) 




22.6 


-1.56(0.03) 


53231.0 


25.36(0.11) 




27.4 


-1.58(0.01) 


53244.4 


25.91(0.14) 




34.3 


-1.58(0.02) 


53284.6 


26.54(0.51) 




54.8 


-1.51(0.02) 




FllOW 






FllOW ^ y 


53208.2 


24.17(0.07) 




15.7 


-1.25(0.03) 


53146.4 


flux:-0.01(0.05) 


-15.8 


(-) 






HST04Rak (z 


= 0.74) 






F606W 






F606W ^ U 


53356.3 


24.15(0.05) 




2.5 


0.362(0.01) 




F775W^ 






F775W ^ B 


53356.3 


22.88(0.02) 




2.5 


-1.11(0.01) 


53373.7 


23.70(0.03) 




12.5 


-1.10(0.01) 


53388.6 


24.56(0.05) 




21.1 


-1.12(0.01) 


53398.1 


25.11(0.07) 




26.5 


-1.12(0.01) 


53404.1 


25.57(0.10) 




30.0 


-1.11(0.01) 


53408.0 


25.64(0.10) 




32.2 


-1.11(0.02) 




F850LP 






F850LP ^ V 
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Table 2 — Continued 



Date" 


Vega Mag 


Epoch(rest) 


K-Corr 


53356.3 


22.96(0.02) 


2.5 


-1.10(0.03) 


53373.7 


23.47(0.02) 


12.5 


-1.04(0.03) 


53378.4 


23.61(0.04) 


15.2 


-1.02(0.05) 


53388.6 


23.92(0.03) 


21.1 


-0.95(0.04) 


53398.1 


24.28(0.05) 


26.5 


-0.90(0.03) 


53404.2 


24.63(0.06) 


30.0 


-0.87(0.02) 


53404.2 


24.42(0.07) 


30.0 


-0.87(0.02) 


53407.9 


24.80(0.08) 


32.2 


-0.86(0.01) 


53313.0 


flux:0.070(0.06) 


-22.3 


(-) 




HST05Zwi {z 


= 0.521) 






F606W 




F606W B 


53415.9 


22.94(0.02) 


5.3 


-0.27(0.03) 


53435.0 


24.13(0.02) 


17.9 


-0.14(0.06) 


53443.1 


24.52(0.10) 


23.2 


-0.04(0.03) 




F775W 




F77bW V 


53403.2 


22.40(0.02) 


-2.9 


-0.86(0.02) 


53415.9 


22.19(0.10) 


5.3 


-0.81(0.04) 


53435.1 


23.04(0.02) 


17.9 


-0.71(0.07) 


53443.2 


23.30(0.02) 


23.3 


-0.61(0.07) 




i^850LP 




i^850LP R 


53403.2 


22.04(0.01) 


-2.9 


-0.84(0.04) 


53415.9 


22.16(0.02) 


5.3 


-0.90(0.04) 


53435.1 


22.95(0.02) 


17.9 


-0.97(0.04) 


53443.2 


22.91(0.02) 


23.3 


-0.90(0.06) 




HST04Haw {z 


= 0.490) 






F606W 




FQQQW B 


53375.7 


25.68(0.04) 


49.4 


-0.16(0.01) 


53427.5 


26.16(0.04) 


84.2 


-0.17(0.01) 


53472.5 


26.82(0.04) 


114.4 


-0.17(0.01) 




F775W 




F77hW V 


53284.0 


23.97(0.03) 


-12.0 


-0.80(0.01) 


53332.9 


23.15(0.02) 


20.7 


-0.71(0.03) 


53375.7 


24.45(0.02) 


49.4 


-0.62(0.02) 


53427.5 


25.08(0.05) 


84.2 


-0.70(0.01) 




i^850LP 




FSbOLP R 


53284.0 


23.81(0.03) 


-12.0 


-0.84(0.05) 


53333.0 


22.90(0.01) 


20.8 


-0.91(0.04) 


53375.7 


23.83(0.05) 


49.4 


-0.82(0.01) 


53427.6 


24.94(0.12) 


84.3 


-0.84(0.01) 


53472.5 


25.82(0.21) 


114.4 


-0.84(0.01) 


53243.9 


flux:0.020(0.08) 


-38.9 


(-) 



Table 2 — Continued 



Date" Vega Mag Epoch(rest) K-Corr 



HST04Kur {z = 0.359) 





F775W 






F775W V 


53355.2 


23.60(0.02) 




-5.3 


-0.70(0.02) 


53366.1 


23.28(0.05) 




2.6 


-0.71(0.02) 


53375.9 


23.74(0.05) 




9.9 


-0.73(0.01) 


53404.0 


24.77(0.04) 




30.5 


-0.85(0.01) 


53415.4 


25.05(0.05) 




38.9 


-0.84(0.01) 


53425.9 


25.03(0.05) 




46.6 


-0.85(0.01) 




FS50LP 






i^850LP R 


53355.3 


23.04(0.02) 




-5.2 


-0.70(0.02) 


53366.1 


22.79(0.02) 




2.6 


-0.68(0.02) 


53375.9 


23.12(0.02) 




9.9 


-0.65(0.01) 


53404.0 


24.09(0.05) 




30.5 


-0.72(0.01) 


53415.5 


24.18(0.06) 




39.0 


-0.70(0.01) 


53425.9 


24.26(0.06) 




46.6 


-0.69(0.02) 






HST04YOW {z 


= 0.46) 






F775W 






F775VF ^ V 


53145.1 


22.73(0.01) 




9.4 


-0.75(0.02) 


53157.0 


23.18(0.01) 




17.6 


-0.72(0.03) 


53169.4 


23.70(0.01) 




26.1 


-0.67(0.02) 


53195.2 


24.42(0.30) 




43.7 


-0.66(0.02) 




F850LP 






F850LP ^ R 


53145.2 


22.53(0.01) 




9.5 


-0.85(0.01) 


53157.1 


23.04(0.01) 




17.7 


-0.86(0.04) 


53169.4 


23.31(0.03) 




26.1 


-0.76(0.02) 


53195.2 


24.06(0.05) 




43.7 


-0.78(0.02) 


53216.3 


24.80(0.09) 




58.2 


-0.80(0.02) 


53221.7 


24.90(0.06) 




61.9 


-0.81(0.01) 


53284.4 


25.49(0.17) 




104.8 


-0.82(0.01) 


53332.7 


25.79(0.20) 




137.9 


-0.82(0.01) 


53097.2 


flux:0.050(0.06) 


-23.3 


(-) 






HST04Man {z 


= 0.854) 






F775W 






F77^W B 


53146.3 


23.38(0.04) 




-1.0 


-1.17(0.03) 


53157.0 


23.55(0.04) 




4.7 


-1.13(0.02) 


53169.4 


24.12(0.05) 




11.4 


-1.11(0.03) 


53195.6 


25.92(0.10) 




25.5 


-0.91(0.01) 




i^850LP 






F850LP ^ V 


53146.3 


23.35(0.05) 




-1.0 


-1.31(0.01) 


53157.1 


23.45(0.02) 




4.8 


-1.28(0.07) 


53169.4 


23.69(0.04) 




11.4 


-1.16(0.07) 
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Table 2 — Continued 



Date" 


Vega Mag 




Epoch(rest) 


K-Corr 


53195.6 


25.04(0.10) 




25.5 


-0.75(0.06) 


53207.7 


25.94(0.21) 




32.0 


-0.68(0.02) 


53098.0 


flux:0.030(0.04) 


-27.0 


(-) 






HST05Spo {z -- 


= 0.839) 






F775W 






F775T^ ^ B 


53427.1 


24.06(0.04) 




14.6 


-1.12(0.04) 


53427.3 


23.95(0.04) 




14.7 


-1.12(0.04) 


53436.6 


24.67(0.05) 




19.8 


-1.04(0.06) 


53436.6 


24.76(0.10) 




19.8 


-1.04(0.06) 


53473.6 


26.30(0.30) 




39.9 


-0.98(0.01) 




F850LP 






F850LP V 


53376.8 


24.78(0.09) 




-12.6 


-1.17(0.04) 


53427.2 


23.66(0.03) 




14.7 


-1.09(0.10) 


53427.4 


23.61(0.04) 




14.8 


-1.09(0.10) 


53436.6 


24.00(0.05) 




19.8 


-0.93(0.10) 


53439.6 


24.17(0.07) 




21.4 


-0.88(0.08) 


53445.5 


24.35(0.06) 




24.6 


-0.80(0.08) 


53454.8 


24.91(0.10) 




29.7 


-0.71(0.03) 


53463.3 


25.11(0.09) 




34.3 


-0.70(0.04) 


53471.7 


25.37(0.12) 




38.9 


-0.73(0.02) 


53473.6 


25.31(0.21) 




39.9 


-0.73(0.03) 






HST04Eag {z 


= 1.02) 






F775W 






F775W U 


53284.9 


23.95(0.05) 




-4.4 


-0.63(0.03) 


53332.1 


26.59(0.60) 




18.9 


-0.84(0.01) 




i^850LP 






F850LP B 


53284.9 


23.76(0.05) 




-4.4 


-1.34(0.02) 


53296.9 


23.60(0.02) 




1.5 


-1.36(0.02) 


53305.6 


23.70(0.03) 




5.8 


-1.37(0.01) 


53316.5 


24.19(0.04) 




11.2 


-1.38(0.02) 


53326.2 


24.62(0.07) 




16.0 


-1.40(0.03) 


53332.2 


24.93(0.11) 




19.0 


-1.42(0.03) 


53334.9 


24.96(0.07) 




20.3 


-1.42(0.03) 




i^llOW^ 






Fiim 


53297.5 


23.57(0.07) 




1.8 


-1.49(0.02) 


53306.3 


23.60(0.07) 




6.1 


-1.45(0.06) 






HST05Gab {z 


= 1.12) 






F775W 






F775VF U 


53474.5 


24.52(0.05) 




-6.0 


-0.50(0.01) 


53485.7 


24.16(0.10) 




-0.7 


-0.50(0.02) 




F850LP 






F850LP B 
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Table 2 — Continued 



Date" 


Vega Mag 


Epoch(rest) 


K-Corr 


53474.5 


23.83(0.04) 


-6.0 


-1.48(0.01) 


53485.7 


23.58(0.03) 


-0.7 


-1.48(0.02) 


53493.9 


23.66(0.03) 


3.0 


-1.48(0.02) 


53504.2 


23.88(0.04) 


7.9 


-1.47(0.02) 


53513.3 


24.24(0.04) 


12.2 


-1.44(0.03) 


53530.8 


25.38(0.13) 


20.5 


-1.38(0.03) 




FllOW 




Fimv V 


53486.4 


23.68(0.07) 


-0.4 


-1.57(0.02) 


53492.3 


23.60(0.07) 


2.3 


-1.57(0.03) 




HST05Str {z 


= 1.01) 






FllbW 




F775W U 


53474.5 


23.93(0.04) 


-0.0 


-0.67(0.08) 


53485.7 


24.25(0.10) 


5.4 


-0.76(0.03) 




F850LP 




F850LP B 


53474.5 


23.67(0.03) 


-0.0 


-1.33(0.04) 


53485.7 


23.90(0.04) 


5.4 


-1.36(0.01) 


53493.9 


24.14(0.05) 


9.5 


-1.36(0.01) 


53504.2 


24.51(0.08) 


14.7 


-1.38(0.03) 


53513.3 


25.08(0.08) 


19.2 


-1.42(0.03) 




FllOW 




FllOW V 


53485.4 


23.85(0.12) 


5.3 


-1.45(0.06) 


53492.2 


23.99(0.12) 


8.7 


-1.41(0.04) 


53426.8 


flux:0.043(0.03) 


-23.8 


(-) 




HST04Sas {z 


= 1.39) 






F775W 




F77?>W U 


148.670 


26.22(0.25) 


-2.4 


0.091(0.14) 




FSbOlp ^ U 






148.700 


24.75(0.07) 


-2.4 


-0.96(0.06) 


156.200 


25.01(0.10) 


0.7 


-0.99(0.07) 


163.600 


25.13(0.08) 


3.8 


-1.04(0.05) 


169.300 


25.74(0.15) 


6.2 


-1.07(0.03) 


176.100 


25.54(0.14) 


9.0 


-1.08(0.03) 


183.500 


25.98(0.20) 


12.1 


-1.10(0.03) 


190.400 


26.29(0.25) 


15.0 


-1.12(0.02) 




FllOW 




FllOW -> B 


164.050 


24.50(0.07) 


4.0 


-1.68(0.05) 


168.060 


24.48(0.07) 


5.6 


-1.69(0.05) 




F160W 




Fimw R 


164.180 


24.06(0.10) 


4.0 


-2.01(0.04) 


169.050 


24.05(0.10) 


6.0 


-2.00(0.03) 


98.6000 


flux:-0.06(0.04) 


-23.3 


(-) 



Table 2 — Continued 



Date" Vega Mag Epoch(rest) K-Corr 



2003aj {z = 1.307) 





F775W 




^^775^" ^ U 


52673.1 


26.62(0.10) 


1.7 


-0.05(0.08) 




J^850LP 




F850LP ^ B 


52673.2 


25.31(0.08) 


1.7 


-1.25(0.03) 


52680.2 


25.53(0.06) 


4.8 


-1.21(0.03) 


52694.5 


26.38(0.16) 


11.0 


-1.16(0.05) 




FUOW 




FllOW V 


52684.2 


25.13(0.11) 


6.5 


-1.49(0.06) 




F160W 




FIQOW R 


52685.1 


24.69(0.12) 


6.9 


-1.83(0.04) 


52627.4 


flux:0.062(0.03) 


-18.0 


(-) 




2002fx {z 


= 1.400) 






F775W 




F775W U 


52495.0 


28.02(1.00) 


-14.2 


0.008(0.07) 


52537.8 


27.12(0.25) 


3.6 


0.300(0.10) 


52580.0 


29.02(1.00) 


21.2 


0.812(0.05) 




F850LP 




F850LP B 


52537.8 


25.21(0.07) 


3.6 


-1.20(0.09) 


52580.5 


27.10(0.27) 


21.4 


-0.82(0.08) 


52490.3 


flux:0.016(0.03) 


-16.1 


(— ) 




2003eq (z 


= 0.84) 






F606W^ 




F606W -> U 


52783.7 


24.58(0.01) 


0.8 


0.707(0.05) 




F775W 




F77^W B 


52783.7 


23.22(0.01) 


0.8 


-1.18(0.03) 


52799.1 


23.67(0.02) 


9.2 


-1.14(0.02) 


52807.3 


24.15(0.04) 


13.6 


-1.12(0.04) 


52819.8 


25.05(0.07) 


20.4 


-1.02(0.06) 


52838.3 


26.22(0.30) 


30.5 


-0.97(0.02) 




i^850LP 




P850LP ^ V 


52783.7 


23.12(0.01) 


0.8 


-1.27(0.02) 


52792.1 


23.28(0.03) 


5.4 


-1.22(0.07) 


52799.1 


23.45(0.02) 


9.2 


-1.16(0.05) 


52807.3 


23.64(0.07) 


13.6 


-1.10(0.09) 


52819.8 


24.30(0.10) 


20.4 


-0.89(0.08) 


52838.3 


25.34(0.12) 


30.5 


-0.69(0.02) 




i^llOW 




FllOW^^ i? 


52792.9 


23.20(0.10) 


5.8 


-1.22(0.01) 


52735.5 


flux:0.087(0.10) 


-25.3 


(-) 



2003es {z = 0.954) 
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Date" 


Vega Mag 


Epoch(rest) 


K-Corr 







FllhW 




F775W ^ U 


52784.5 


24.12(0.03) 


8.0 


-0.87(0.05) 




F850LP 




F850LP B 


52784.5 


23.73(0.02) 


8.0 


-1.33(0.03) 


52792.4 


24.09(0.06) 


12.0 


-1.36(0.06) 


52801.3 


24.68(0.06) 


16.6 


-1.42(0.07) 


52807.9 


24.98(0.09) 


19.9 


-1.49(0.06) 


52821.0 


25.63(0.14) 


26.6 


-1.56(0.01) 


52838.1 


25.96(0.24) 


35.4 


-1.58(0.02) 




FllQW 




FllQW^V 


52792.8 


24.20(0.08) 


12.2 


-1.28(0.05) 


52734.6 


flux:0.005(0.04) 


-17.5 


(-) 




2003az {z = 


1.265) 






F775W 




F77hW U 


52690.9 


25.10(0.05) 


3.0 


-0.15(0.06) 


52701.2 


25.46(0.05) 


7.5 


-0.09(0.05) 




F850LP 




F850LP B 


52690.9 


24.36(0.04) 


3.0 


-1.40(0.03) 


52701.2 


24.48(0.04) 


7.5 


-1.37(0.02) 


52709.1 


24.64(0.05) 


11.0 


-1.35(0.03) 


52716.9 


25.09(0.06) 


14.4 


-1.32(0.05) 


52726.5 


25.52(0.08) 


18.7 


-1.25(0.08) 


52733.2 


25.75(0.09) 


21.6 


-1.17(0.07) 




FllOW 




FllQW 


52703.6 


24.10(0.06) 


8.6 


-1.51(0.06) 


52710.6 


24.25(0.06) 


11.7 


-1.46(0.07) 


52642.2 


flux:0.054(0.03) 


-18.4 


{-) 




2002kc {z = 


0.216) 






F6mw 




F&Q&W V 


52629.6 


22.30(0.01) 


-7.4 


-0.23(0.05) 


52672.3 


23.28(0.01) 


27.6 


0.222(0.02) 




F775W 




F77bW R 


52629.6 


21.81(0.02) 


-7.4 


-0.43(0.01) 


52642.5 


21.35(0.10) 


3.1 


-0.47(0.04) 


52672.3 


22.11(0.01) 


27.6 


-0.43(0.01) 




F850LP 




F850LP ^ / 


52629.7 


21.70(0.01) 


-7.4 


-0.45(0.01) 


52672.3 


21.87(0.01) 


27.6 


-0.19(0.07) 




2003eb {z = 


= 0.90) 






F606W 




FQmw U 


52734.6 


24.26(0.02) 


-1.1 


0.824(0.07) 
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Date" 


Vega Mag 


Epoch(rest) 


K-Corr 


52783.5 


27.07(0.15) 


24.5 


1.233(0.02) 




F775W 




F775W B 


52734.6 


23.05(0.02) 


-1.1 


-1.14(0.05) 


52745.6 


23.15(0.02) 


4.5 


-1.10(0.03) 


52783.5 


25.22(0.06) 


24.5 


-0.83(0.02) 


52799.1 


25.92(0.10) 


32.7 


-0.80(0.03) 




^^850LP 




F8FMLP V 


52734.7 


22.83(0.01) 


-1.1 


-1.32(0.01) 


52745.7 


22.84(0.01) 


4.6 


-1.27(0.09) 


52751.2 


22.98(0.01) 


7.5 


-1.20(0.06) 


52763.6 


23.47(0.01) 


14.0 


-1.09(0.14) 


52774.0 


24.02(0.04) 


19.5 


-0.92(0.13) 


52783.6 


24.39(0.04) 


24.5 


-0.75(0.10) 


52792.1 


24.77(0.11) 


29.0 


-0.60(0.03) 


52792.4 


24.63(0.08) 


29.2 


-0.60(0.03) 


52799.1 


25.06(0.15) 


32.7 


-0.55(0.02) 


52801.3 


24.84(0.06) 


33.9 


-0.56(0.04) 


52807.9 


25.03(0.10) 


37.3 


-0.61(0.02) 


52821.0 


25.40(0.09) 


44.2 


-0.63(0.02) 


52838.1 


25.36(0.10) 


53.2 


-0.69(0.05) 


52692.5 


flux:0.058(0.05) 


-23.3 


(-) 




2003XX {z 


= 0.935) 






F606W 




F606W U 


52733.6 


25.17(0.05) 


5.7 


1.042(0.09) 




F775W 




F775W B 


52733.6 


23.57(0.05) 


5.7 


-1.08(0.02) 


52745.6 


24.12(0.05) 


11.9 


-1.05(0.04) 


52783.5 


26.40(0.25) 


31.5 


-0.77(0.02) 




F850LP 




F850LP ^ V 


52692.5 


26.20(0.13) 


-15.4 


-1.12(0.06) 


52733.7 


23.26(0.01) 


5.7 


-1.27(0.10) 


52745.7 


23.63(0.01) 


11.9 


-1.14(0.10) 


52751.2 


23.88(0.02) 


14.8 


-1.07(0.15) 


52763.6 


24.56(0.04) 


21.2 


-0.77(0.11) 


52774.0 


24.87(0.05) 


26.6 


-0.58(0.05) 


52783.6 


25.20(0.06) 


31.5 


-0.48(0.02) 


52792.1 


25.36(0.13) 


35.9 


-0.50(0.05) 


52799.1 


26.04(0.20) 


39.5 


-0.53(0.04) 


52807.3 


26.45(0.19) 


43.8 


-0.56(0.02) 




20021ir (- : 


= 0.32()) 






F606W 




FQmw u 
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Table 2 — Continued 



Date" 


Vega Mag 




Epoch(rest) 


K-Corr 


52579.6 


24.04(0.03) 




-6.4 


0.103(0.08) 


52629.5 


25.82(0.15) 




26.2 


-0.36(0.03) 


52674.1 


27.43(0.30) 




55.5 


-0.35(0.07) 




F775W 






F775W B 


52579.6 


23.56(0.03) 




-6.4 


-0.84(0.01) 


52589.7 


23.20(0.12) 




0.1 


-0.87(0.05) 


52590.0 


23.28(0.11) 




0.3 


-0.87(0.05) 


52596.6 


23.44(0.14) 




4.7 


-0.93(0.09) 


52614.5 


23.91(0.08) 




16.4 


-1.27(0.18) 


52629.5 


24.42(0.05) 




26.2 


-1.66(0.03) 


52674.1 


25.67(0.10) 




55.5 


-1.58(0.02) 




F850LP 






F850LP V 


52579.7 


23.38(0.02) 




-6.3 


-0.96(0.05) 


52629.5 


23.76(0.03) 




26.2 


-1.18(0.04) 


52674.1 


25.04(0.10) 




55.5 


-1.14(0.02) 






2003bd {z = 


= 0.67) 






F606W 






Fmew u 


52691.9 


24.57(0.04) 




8.1 


0.162(0.02) 


52735.4 


27.50(0.12) 




34.1 


0.170(0.02) 




F775W 






F77bW B 


52691.9 


23.41(0.03) 




8.1 


-1.05(0.02) 


52735.4 


25.77(0.10) 




34.1 


-1.24(0.02) 


52745.6 


26.05(0.10) 




40.2 


-1.25(0.01) 




F850LP 






F850LP ^ V 


52692.0 


23.18(0.01) 




8.1 


-1.07(0.01) 


52735.5 


24.82(0.08) 




34.2 


-1.02(0.01) 


52745.7 


25.14(0.08) 




40.3 


-1.02(0.01) 


52751.2 


25.21(0.09) 




43.6 


-1.02(0.01) 


52763.6 


25.47(0.14) 




51.0 


-1.03(0.01) 


52774.0 


25.56(0.19) 




57.2 


-1.04(0.01) 


52792.1 


25.96(0.21) 




68.1 


-1.04(0.01) 


52807.3 


25.85(0.26) 




77.2 


-1.04(0.01) 


52819.8 


25.88(0.27) 




84.6 


-1.04(0.01) 


52838.3 


26.57(0.31) 




95.7 


-1.04(0.01) 


52642.6 


flux:0.057(0.03) 




-21.4 


(-) 






2002kd {z = 


0.735) 






F606W 






Fm&w u 


52629.4 


24.96(0.03) 




-8.9 


0.335(0.01) 


52673.5 


25.81(0.05) 




16.4 


0.335(0.03) 




F775W 






F775W B 


52629.4 


23.78(0.02) 




-8.9 


-1.06(0.04) 
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Date" 


Vega Mag 




Epoch(rest) 


K-Corr 


52644.5 


22.96(0.17) 




-0.2 


-1.10(0.03) 


52645.5 


22.76(0.10) 




0.3 


-1.10(0.03) 


52673.5 


24.18(0.05) 




16.4 


-1.11(0.02) 




F850LP 






F850LP V 


52629.5 


23.76(0.02) 




-8.8 


-1.10(0.01) 


52639.4 


23.09(0.01) 




-3.1 


-1.12(0.03) 


52673.7 


23.68(0.02) 




16.5 


-1.00(0.05) 






2003be {z ■ 


= 0.640) 






F606W 






P606PF -> U 


52692.0 


24.47(0.03) 




12.1 


0.096(0.03) 


52732.5 


27.26(0.20) 




36.8 


0.086(0.03) 




FllbW 






FllbW B 


52641.3 


28.02(0.80) 




-18.7 


-0.68(0.12) 


52692.0 


23.45(0.03) 




12.1 


-1.04(0.06) 


52732.5 


25.50(0.10) 




36.8 


-1.29(0.02) 


52784.3 


26.09(0.10) 




68.4 


-1.22(0.02) 




F850i:P 






P850i:P V 


52692.1 


23.04(0.01) 




12.2 


-1.07(0.01) 


52732.5 


24.45(0.05) 




36.8 


-1.08(0.01) 


52784.4 


25.31(0.13) 




68.4 


-1.07(0.01) 


52641.3 


flux:0.087(0.10) 




-18.7 


(-) 






2003dy [z 


= 1.34) 






F850LP 






F850i:P ^ u 


52733.7 


24.45(0.04) 




-0.1 


-0.96(0.09) 


52745.7 


24.67(0.06) 




4.9 


-1.05(0.05) 


52751.2 


24.77(0.08) 




7.3 


-1.09(0.04) 


52763.6 


25.29(0.11) 




12.6 


-1.16(0.04) 


52774.0 


25.57(0.15) 




17.0 


-1.21(0.01) 


52783.6 


26.62(0.45) 




21.1 


-1.22(0.02) 


52801.3 


26.65(0.26) 




28.7 


-1.27(0.03) 




FllQW 






FllQW B 


52751.6 


24.28(0.07) 




7.4 


-1.73(0.05) 


52754.6 


24.38(0.08) 




8.7 


-1.74(0.05) 




F160W 






Fimw R 


52751.7 


23.81(0.08) 




7.5 


-1.89(0.01) 


52692.4 


flux:0.037(0.04) 




-17.8 


(-) 






2002ki {z 


= 1.14) 






F775W^ 






FllbW -> U 


52600.8 


24.79(0.10) 




-0.6 


-0.46(0.01) 


52643.6 


26.82(0.17) 




19.3 


-0.46(0.01) 




F850LP 






P850LP B 
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Date" 


Vega Mag 




Epoch(rest) 


K-Corr 


52600.8 


23.89(0.04) 




-0.6 


-1.48(0.02) 


52643.6 


25.80(0.10) 




19.3 


-1.36(0.05) 


52652.2 


25.99(0.10) 




23.3 


-1.32(0.02) 


52663.7 


26.89(0.20) 




28.7 


-1.32(0.02) 




F160W 






F160W I 


52652.4 


24.68(0.25) 




23.4 


-1.65(0.03) 


52664.5 


24.87(0.25) 




29.1 


-1.55(0.07) 






2003ak {z = 


1.551) 






F850LP 






FSbOLP U 


52673.2 


25.68(0.14) 




5.0 


-0.68(0.03) 


52680.2 


25.95(0.10) 




7.7 


-0.67(0.03) 


52694.5 


26.62(0.21) 




13.3 


-0.63(0.03) 




FUOW 






FllOW ^ B 


52681.1 


25.03(0.10) 




8.1 


-1.71(0.03) 


52693.0 


25.03(0.10) 




12.8 


-1.74(0.05) 


52708.5 


25.38(0.15) 




18.8 


-1.82(0.07) 


52715.4 


25.80(0.15) 




21.5 


-1.87(0.04) 




F160W^ 






F160W^ V 


52681.3 


24.01(0.05) 




8.2 


-2.29(0.03) 


52693.4 


24.21(0.07) 




12.9 


-2.25(0.03) 


52701.6 


24.61(0.08) 




16.1 


-2.26(0.11) 


52627.4 


flux:0.050(0.03) 




-12.9 


(-) 






2002hp {z = 


1.305) 






F775W 






F775W U 


52579.5 


25.52(0.10) 




3.3 


-0.03(0.07) 




^^850LP 






F850LP B 


52537.1 


26.83(0.39) 




-15.0 


-1.47(0.06) 


52579.5 


24.32(0.04) 




3.3 


-1.27(0.05) 


52589.1 


24.89(0.06) 




7.4 


-1.23(0.03) 


52595.4 


25.06(0.08) 




10.1 


-1.22(0.03) 


52603.8 


25.67(0.08) 




13.8 


-1.19(0.05) 


52613.8 


26.09(0.14) 




18.1 


-1.10(0.09) 


52629.3 


26.60(0.17) 




24.9 


-1.00(0.03) 


52639.4 


27.30(0.35) 




29.2 


-0.99(0.04) 




FllOW 






FllQW V 


52589.2 


24.28(0.07) 




7.5 


-1.52(0.06) 


52595.5 


24.33(0.07) 




10.2 


-1.46(0.06) 






2002fw {z = 


: 1.30) 






F775W 






F775W U 


52536.8 


25.37(0.05) 




-5.8 


-0.22(0.04) 


52548.3 


25.01(0.08) 




-0.8 


-0.14(0.07) 
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Table 2 — Continued 



Date" 


Vega Mag 




Epoch(rest) 


K-Corr 


52578.4 


26.45(0.10) 




12.2 


0.111(0.08) 




F850LP 






F850i:P ^ B 


52536.9 


24.49(0.05) 




-5.8 


-1.40(0.02) 


52548.3 


24.17(0.04) 




-0.8 


-1.40(0.07) 


52552.5 


24.44(0.06) 




0.9 


-1.40(0.07) 


52557.4 


24.36(0.05) 




3.1 


-1.36(0.05) 


52567.9 


24.48(0.06) 




7.6 


-1.33(0.02) 


52577.5 


24.85(0.07) 




11.8 


-1.31(0.04) 


52578.5 


25.07(0.07) 




12.2 


-1.30(0.05) 


52595.5 


25.88(0.17) 




19.6 


-1.17(0.08) 


52603.8 


26.68(0.25) 




23.2 


-1.07(0.03) 




FllOW 






FllOW V 


52549.5 


23.96(0.08) 




-0.3 


-1.67(0.04) 


52557.6 


24.11(0.08) 




3.1 


-1.64(0.07) 




F160W 






F160W R 


52549.8 


23.83(0.08) 




-0.2 


-1.92(0.09) 


52557.8 


23.76(0.09) 




3.2 


-1.85(0.07) 






2002dc {z 


= 0.475) 






F775W 






F775T^ B 


52405.4 


22.41(0.02) 




-0.0 


-0.78(0.05) 


52438.0 


23.57(0.03) 




22.0 


-1.66(0.13) 


52449.3 


24.07(0.05) 




29.7 


-1.78(0.02) 




F850LP 






F850LP V 


52415.2 


22.38(0.02) 




6.6 


-0.89(0.06) 


52438.1 


23.23(0.03) 




22.1 


-1.02(0.08) 


52449.4 


23.51(0.06) 




29.7 


-1.14(0.04) 


52455.2 


23.85(0.05) 




33.7 


-1.15(0.05) 


52600.0 


27.54(0.65) 




131.8 


-0.94(0.01) 






2002dd {z 


= 0.95) 






F775W 






F775W -> U 


52405.4 


23.66(0.02) 




0.8 


-0.66(0.12) 


52438.0 


24.99(0.06) 




17.5 


-0.98(0.01) 


52449.3 


25.51(0.09) 




23.3 


-1.01(0.03) 




F850LP 






F850LP B 


52405.4 


23.24(0.05) 




0.8 


-1.27(0.05) 


52415.2 


23.45(0.02) 




5.8 


-1.31(0.02) 


52438.1 


24.34(0.08) 




17.6 


-1.44(0.08) 


52449.4 


24.85(0.12) 




23.4 


-1.56(0.02) 


52455.2 


24.96(0.09) 




26.3 


-1.58(0.01) 




FllOW 






FllOW 


52438.0 


23.98(0.05) 




17.5 


-1.27(0.03) 
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3. Distance Fitting and Cosmology Constraints 



Distance estimates from SN la light curves are derived from the luminosity distance, 



where £ and T are the intrinsic luminosity and the absorption-free flux of the SN 
within a given passband, respectively. Equivalently, logarithmic measures of the flux 
(apparent magnitude, m), luminosity (absolute magnitude, Af), and colors (to quantify 
the selective absorption) were used to derive extinction-corrected distance moduli, 
Ho = m — M = blogd^ + 25 {d^ in units of megaparsecs) . In this context, the luminosity is 
a "nuisance parameter" whose fiducial value vanishes from a reconstruction of the expansion 
history (which makes use of differences in distance with redshift). We have adopted the 
MLCS2k2 method (Riess, Press, & Kirshner 1996; Jha 2002) used by R04 and the data 
described in §2 to derive accurate and individual relative distance moduli for the sets of SNe 
described in §3 and given in Table 6^^. As in Riess, Press, & Kirshner (1996), we require 
for robust fits of light curves that SN photometry commence no less than 10 days after 
maximum, although in practice degeneracies in light curve fits for late-commencing SN la 
photometry are also alleviated by flux limits from the preceding epoch (typically 3 weeks 
earlier than the discovery point in the restframe). 



^® Although a new version of the MLCS2k2 algorithm is currently under development with an expanded 
training set (Jha, Riess, & Kirshner 2006), we utilized the same version used by R04 to maintain, as far as 
possible, a consistent and tested process to measure all available SN la data. However, tests comparing 
the R04 and Jha et al. (2006) versions of the algorithm show very good agreement when applied to 
high quality data (Lampeitl 2006). In addition we have found and corrected a minor numerical error 
in our calculation of ii'-corrections which affect a few SNe la at the few percent level. In the future, 
it will be valuable to reanalyze all available data consistently with continually improving algorithms and 
we will attempt to provide updated distance esimates using updated versions of data and algorithms at 
http://braeburn.pha.jhu.edu/~ariess/R06 or upon request to ariess@stsci.edu. 
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Table 2 — Continued 



Date' 



Vega Mag 



Epoch(rest) 



K-Corr 



^JD-2,400,000. 
Uncertainties in magnitudes are listed in parentheses. 
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Table 3. 


Spectroscopic Data 




SN 


MJD(age) 


Instrument 


z 





HST04Pat 


53205.1(+3) 


HST KQ.^ 


0.97"'^ 


HST04Cay 


53254.9 


EST ACS 




HST04Mcg 


53275.6(+10) 


ifSTACS 


1.357'='i 


HST05Fer 


53486.8(+10) 


ESTKQ.'^ 


1.02°'i 


HST05Koe 


53486.8(+15) 


HST ACS 


1.23"'i 


HST05Dic 






0.638'''3 


HST04Gre 


53275.6(+6) 


EST ACS 


1.14a,i 


HST04Oinb 


53373.7(+3) 


EST ACS 


0.975a,d,i 


HSTOSRed 


53436.6 


EST ACS 


i.isg''^^ 


HSTOSLan 


53436.6(+2) 


EST ACS 


1.235'='^ 


HST04Tha 


53207.7(+15) 


EST ACS 


0_954a,d,i 


HST04Rak 


53373.7(+12) 


EST ACS 


0.739^'^ 


HST05Zwi 


53415.9(+5) 


EST ACS 




HST04Haw 






0.490'^'^ 


HST04Kur 


53366. l(+3) 


EST ACS 


0.359^'^ 


HST04YOW 


53157.0(+17) 


EST ACS 


Q 457a,d,i 


HST04Man 


53157.0(+5) 


EST ACS 


0.854"''''! 


HST05Spo 


53436.6(+20) 


EST ACS 


0.839"''''i 


HST04Eag 


53296.9(+2) 


EST ACS 


1.019'='^ 


HST05Gab 


53485.7(-l) 


EST ACS 


1.12°'i 


HST05Stro 


53485. 7(+5) 


EST ACS 


2 027°''''i 


HST04Sas 


53156.2(+1) 


EST ACS 


1.39"'^ 



''Prom cross-correlation with broad SN features. 

''Prom narrow features in our own additional Keck host-galaxy spectrum. 
'^Prom both (a) and (b). 
'^From Team Keck. 

"'Classified as SN la with high confidence from spectrum. 
^Photometric properties indicate likely (but not certain) SN la. 
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In Figure 6 we show the Hubble diagram of distance moduh and redshifts for all of the 
if^T-discovered SNe la in the Gold and Silver sets from our program. The new SNe la span 
a wide range of redshift (0.21 < z < 1.55), but their most valuable contribution to the SN la 
Hubble diagram remains in the highest-redshift region where they now well delineate the 
range at 2; > 1 with 23 SNe la, 13 new objects since R04. This territory remains uniquely 
accessible to HST, which has discovered the dozen highest-redshift SNe la known, and its 
exploration is the focus of the rest of this paper. 

In the inset to Figure 6 we show the residual Hubble diagram (from an empty Universe) 
with the Gold data uniformly binned. Here and elsewhere, we will utilize uniform, unbiased 
binning achieved with a fixed value of nAz, where Az is the bin width in redshift and n is 
the number of SNe in the bin.^^ In Figure 6 we use nAz — 6 which yields seven bins for our 
sample. Although binning is for illustrative purposes in the Hubble diagram, there are some 
specific advantages of binning such as the removal of lensing-induced asymmetrical residuals 
by flux averaging (Wang 2005) and the ease of accounting for systematic uncertainties 
introduced by zeropoint errors in sets of photometric passbands used at similar redshifts. 

The distance-redshift relation of SNe la is one of few powerful tools available in 
observational cosmology. A number of different hypotheses and models can be tested with 
it, including kinematic descriptions of the expansion history of the Universe, the existence 
of mass-energy terms on the right-hand side of the Friedman equation, and the presence of 
astrophysical sources of contamination. Testing all interesting hypotheses is well beyond 
the scope of this paper and is best left for future work. Instead, we now undertake a few 
narrowly posed investigations. 

For the following analyses we limit the low-redshift boundary of our sample to cz > 7000 
km s^^ to avoid the influence of a possible local, "Hubble Bubble" (Jha 2002, Zehavi et al. 
1998). The Gold sample above this velocity boundary consists of 182 supernovae. 

3.1. Kinematics and Model-Independent Observables 

The distance-redshift relation of SNe la contains a cumulative record of changes in 
the cosmic expansion rate, H{z), over the last 10 Gyr, described in flat space as given in 
equation 3. Although the interesting kinematic information, H{z), appears to be degraded 
by the physical integration across its temporal variations, good sampling of SNe la over 
a wide range of redshifts allows us to recover its value at discrete, uncorrelated epochs. 



^The last bin ends abruptly with the highest-redshift SN; thus, its nAz < value is smaller than the rest. 
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independent of the cosmological model. Such information may be more general and of 
longer-lasting value than constraints on any single, specific model of dark energy. 

Following Wang & Tegmark (2005; see also Daly & Djorgovski 2004), we transform the 
Gold sample of luminosity distances to comoving coordinate distances, r{z), as 

2997.9(1 + .) ^»"''" 

We assume spatial flatness (as motivated by theoretical considerations, i.e., that most 
inflation models predict Qk < 10~^ or by similar resolutions to the "flatness problem") for 
ease of calculations, but the following approach can be generalized to allow for nontrivial 
spatial curvature. After sorting the SNe la by redshift, we define the quantity 

ri+i — Ti 

Xi = , (6) 

Zi+l — Zi 

where the mean value of Xi gives an unbiased estimate of the inverse of H{z) at the redshift, 
Zi. As in Wang & Tegmark (2005), wc flux average the data first {n/\z ^ 1) to remove 
possible Icnsing bias. Wc then calculate the minimum-variance values of H{z) in 3, 4, or 
5 even-sized bins across the sample, with n^z chosen to be 40, 20, or 15, respectively, to 
achieve the desired number of bins. 

In the upper panel of Figure 7 we show sets of 3, 4, or 5 samplings of H{z) versus 
redshift from the Gold sample. As seen, H[z) remains well constrained until z ~ 1.3, 
beyond which the SN sample is too sparse to usefully determine H{z). For comparison, 
we show the dynamical model of H{z) derived from H(z)'^ = ifQ(f2M(l + z)^ + VIa) with 
"concordance" values of Qm — 0.29 and Qa = 0.71. 

In the lower panel of Figure 7 we show the kinematic quantity a = H{z)/{1 + z) versus 
redshift. In the uncorrelated d versus redshift space, it is very easy to evaluate the sign of 
the change in expansion rate independent of the cosmological model. For comparison we 
show three simple kinematic models: purely accelerating, decelerating, and coasting, with 
q{z) = {-a/a)/H'^{z) = d{H-^{z))/dt - 1 = 0.5, -0.5, and 0.0, respectively. We also show 
a model with recent acceleration (go = —0.6) and previous deceleration dq/dz — 1.2, where 
q{z) — qo + zdq/dz, which is a good fit to the data. 

In Figure 8 we demonstrate the improvement in the measure oi H{z) at z > 1 realized 
from the addition of the new SNe la, presented here, to the sample from R04: we have 
reduced the uncertainty of H{z) at 2; > 1 from just over 50% to just under 20%.^^ 



^^Monte Carlo simulations of the determination of uncorrelated components of H{z) show that the increase 
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We also repeated the analysis of R04 in which the deceleration parameter, 

q{z) = {—a/a)/H'^(z) = d{H~^{z))/dt — 1, is parameterized by q{z) = qo + zdq/dz and 
determined from the data and Eq. (5). As in R04, we find that the Gold set strongly favors 
a Universe with recent acceleration (go < 0) and previous deceleration [dq/dz > 0) with 
99.96% and 99.999% likelihood (summed within this quadrant), respectively Summing the 
probability density in the qo vs. dq/dz plane along lines of constant transition redshift, 
= —%/ (dq/dz), yields the transition redshift of Zt = 0.43 ± 0.07. 

However, as shown by Shapiro & Turner (2005), different parameterizations of q{z) 
can lead to different rcdshifts for the transition and to different confidence levels for the 
epochs of acceleration and deceleration (though all appear sufficiently high to yield a robust 
conclusion similar to our own; see also Daly & Djorgovski 2004). Thus, our uncertainty 
is only statistical within a linear form for q{z), and is not useful for comparing to the 
transition expected for a cosmological constant (zt ~ 0.7) due to the different functional 
forms used for evaluating the redshift at which d — 0. 

3.2. Alternatives to Dark Energy 

After the detection of the apparent acceleration of cosmic expansion (and dark energy) 
by Riess et al. (1998) and Perlmutter et al. (1999), alternative hypotheses for the apparent 
faintness of high-redshift SNe la were posed. These included extragalactic gray dust with 
negligible tell-tale reddening or additional dispersion (Aguirre 1999a,b; Rana 1979, 1980), 
and pure luminosity evolution (Drell, Loredo, & Wasserman 2000). 

As reported elsewhere, there is considerable evidence against these possibilities and 
little evidence in favor of either of them (Riess et al. 1998; Perlmutter et al. 1999; Riess 
et al. 2001; Leibundgut 2001; Sullivan ct al. 2003; Knop et al. 2003; Filippenko 2004, 
2005). However, it is important to remain vigilant for such possibilities. The redshift range 
z > 1, where deceleration dominates over acceleration for the simplest cosmological models, 
is crucial for breaking degeneracies between astrophysical effects and cosmological effects 
(e.g., Riess ct al. 2001). 

In R04 we found that the first significant sample of SNe la at 2; > 1 from i/S'T rejected 
with high confidence the simplest model of gray dust by Goobar, Bergstrom, & Mortsell 
(2002), in which a smooth background of dust is present (presumably ejected from galaxies) 



in precision proceeds as approximately n^^^, significantly faster than n^/^, where n is the number of SNe 
due to the rate of increase in unique pairs of SNe. 
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at a redshift greater than the SN sample (i.e., z > 2) and diluted as the Universe expands. 
This model and its opacity was invented to match the 1998 evidence for dimming of 
supernovae at z 0.5 without invoking dark energy in a universe with Q„i = 1- This model is 
shown in the inset of Figure 6. The present Gold sample (at the best fitting value of Hq) 
rejects this model at even higher confidence (A^^ = 194, i.e., 14 a, see Table 4), beyond a 
level worthy of further consideration. 

However, a more pernicious kind of dust was also suggested by Goobar, Bergstrom, & 
Mortscll (2002), a "replenishing dust" in which a constant density of gray dust is continually 
replenished at just the rate it is diluted by the expanding Universe. This latter dust is 
virtually indistinguishable from an Q/^ model (see Table 4) via the magnitude-redshift 
relation because the dimming is directly proportional to distance traveled and thus 
mathematically similar to the effects of a cosmological constant. Dust of this sort with the 
required opacity, replenishing rate, and ejection velocity from galaxies (> 1000 km s"^ for 
it to fill space uniformly without adding detectable dispersion) may always be virtually 
undetectable in the Hubble diagram of SNe la, but its degree of fine tuning makes it 
unattractive as a simple alternative to a cosmological constant. 

More recently, Ostman & Mortsell (2005) used 11,694 quasars from the Sloan Digital 
Sky Survey (SDSS), at 0.1 < ^ < 2, to hmit the dimming of SNe la at ^ = 1 by Milky Way 
dust to less than 0.03 mag, and by very gray dust (-Rv' = 12) to less than 0.1 mag, for either 
the high-z dust or replenishing dust models. Petric et al. (2006) has ruled out the presence 
of intergalactic gray dust at a level of Qdust < 10~^ by the lack of an X-ray scattering halo 
around a quasar a.t z = 4.3, limiting dimming due to micron-sized grains to less than a few 
percent (less than a percent to 2; ~ 0.5) and this dust scenario has also been limited by the 
resolution of the far infrared background (Aguirre & Haiman 2000) 

We also find, as in R04 but with even greater confidence, that the data are inconsistent 
(Ax^ = 116) with the simplest SN luminosity evolution (proportional to redshift) in lieu of 
a cosmological constant; see inset of Figure 6 and Table 4. 

Less direct but significant failures to detect SN evolution come from detailed 
comparisons of the composite state of the abundance, temperature, and outward velocity 
of the SN photosphere as recorded through the spectra (Sullivan 2005, Blondin et al. 2006, 
Hook et al. 2005, Balland et al. 2006). Greater leverage for this test can be gained by 
extending the redshift range of the observed SED using the present sample. In Figure 16 
we show the average, composite spectrum of the 13 best-observed SNe la at 2; > 1. For the 
spectra listed in Table 3 (and from R04), we transformed to the rest frame and calculated 
a 3(T-clipped average at each wavelength point. Shown are the mean and dispersion about 
the mean. For comparison, we used the template spectra from Nugent et al. (2002) (with 
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colors matched to the MLCS empirical model) averaged over the week following maximum 
light. As shown, the high-redshift composite, with a mean redshift of 1.2, bears a striking 
resemblance to the model of a post-maximum SN la. Two specific indications of SNe la 
are seen in the Si feature at 4130 A and the width of the Ca II and Si II blended feature 
at 3750 A which indicates the presence of both ions. We put less reliance on the region 
redward of 4300 A corresponding to rcdward of observed 9500 A because 1) the quantum 
efficiency of the detector is rapidly decaying and thus harder to calibrate and 2) modest 
color differences may exist between the spectral template and the mean of these objects 
(due to intrinsic variation or reddening). 

We find no evidence for a difference in the mean SED of SNe la across 10 Gyr of 
look-back time or 1.2 units of redshift. 

SNe la hosted by red, carly-typc hosts are expected to be particularly robust against 
absorption by dust (where the dust content is low) and may be expected to arise from 
the earliest forming progenitors (as compared to late- type hosts). Thus differences in 
distances measured to SNe la hosted by early and late-type hosts form an important probe 
of systematics (Sullivan et al. 2003). We find that the 6 SNe la we discovered at z > 1 with 
ACS residing in red, elliptical hosts (HSTOSLan, 2002hp, HST04Sas, 2003es, HST04Tha, 
and 2003XX) all have low measured extinction {Ay < 0.25 mag) and the same dispersion 
and mean (closer for their redshift by 0.07 ±0.10 mag) as the full sample, indicating that 
at 2; > 1 no differences are apparent for distances measured in early type hosts. 



3.3. Dark Energy 

Strong evidence suggests that high-redshift SNe la provide accurate distance 
measurements and that the source of the apparent acceleration they reveal lies in the 
negative pressure of a "dark energy" component. Proceding from this conclusion, our 
hard-earned sample of SNe la at 2; > 1.0 can provide unique constraints on its properties. 
Strong motivation for this investigation comes from thorough studies of high-redshift 
and low-redshift SNe la, yielding a consensus that there is no evidence for evolution 
or intergalactic gray dust at or below the current statistical constraints on the average 
high-redshift apparent brightness of SN la (see Filippenko 2004, 2005 for recent reviews). 
We summarize the key findings here. (1) Empirically, analyses of SN la distances (after 
accounting for the light curve shape-luminosity relation) versus host stellar age, morphology 
or dust content (Riess et al. 1998, Sullivan et al. 2003, Jha, Riess, & Kirshner 2006), or 
metallicity or star formation rate (Gallagher et al. 2005, Riess et al. 1999) indicate that 
SN la distances are relatively indifferent to the evolution of the Universe. (2) Detailed 
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examinations of the distance-independent properties of SNe la (including the far-UV flux, 
e.g., as presented in the last section) provide strong evidence for uniformity across rcdshift 
and no indication (thus far) of rcdshift-dependent differences (e.g., Sullivan et al. 2005; 
Howell et al. 2005; Blondin et al. 2006). (3) SNe la are uniquely qualified as standard 
candles because a well understood, physical limit (the Chandrasekhar limit) provides the 
source of their homogeneity. Based on these studies, we adopt a limit on redshift-dependent 
systematics is to be 5% per /S.z = 1 at 2; > 0.1 and make quantitative use of this in §4.1. 

Many have studied the constraint placed by the rcdshift-magnitude relation of SNe la 
on the parameter combination VLm-w^ where w (assumed to be constant) is the dark energy 
equation-of-state parameter. There are few models for dark energy that predict an equation 
of state that is constant, different from the cosmological constant, and not already ruled 
out by the data. On the other hand, a prominent class of models does exist whose defining 
feature is a time-dependent dark energy (i.e., quintessence). While the rejection oiw — —l 
for an assumed constant value of w would invalidate a cosmological constant, it is also 
possible that apparent consistency with = — 1 in such an analysis would incorrectly 
imply a cosmological constant. For example, if w{z) is rising, declining, or even sinusoidal, 
a measured derivative could be inconsistent with zero while the average value remains 
near —1. Therefore, when using w^z) to discriminate between dark-energy models, it is 
important to allow for time-varying behavior, or else valuable information may be lost. 
Here, we seek to constrain the value of w{z > 1) and bound its derivative across the range 
0.2 < 2; < 1.3. This is unique information afforded by the if^T-discovered SN la sample. 

Unfortunately, present dynamical dark-energy models in the literature (see Szydlowski, 
Kurek, & Krawiec 2006 for a review) do not suggest a universal or fundamental parametric 
form for w{z). Instead, most models contain embedded, free-form functions (e.g., the shape 
of a scalar potential). Thus, we proceed with simple, empirical studies of variations in the 
equation-of-state parameter of dark energy. 

The luminosity distance to SNe la, through the solution to the Preidman equation, is 
sensitive to w{z): 



Here, E{z) = {nm{l + z'f + {I - x exp[+3 C^^-^'' d\n{l + z'){l + w{z'))]Y/\ where a 



is the dimensionless matter density Sti pm / {"^Hq) , Hq is the Hubble constant (the present 
value of the Hubble parameter), and z is the redshift of any SN la. 

We can determine the likehhood for parameters of w{z) from a statistic. 




(7) 



x2(i/o,^^M,Wi) = J2 
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where is the dispersion in supernova redshift (transformed to units of distance moduli) 
due to peculiar velocities a^^ . is the uncertainty in the individual distance moduli, and 
Wi is a set of dark-energy parameters describing wi^z). Due to the extreme redshift of 
our distant sample and the abundance of objects in the nearby sample, our analysis is 
insensitive to the value we assume for within its likely range of 200 km s^^ < o-^, < 
500 km s~^. For our analysis we adopt = 400 km s~^ For high-redshift SNe la whose 
redshifts were determined from the broad features in the SN spectrum, we assume an 
uncertainty of — 2500 km in quadrature to the peculiar velocity term. Marginalizing 
our likelihood functions over the nuisance parameter, Hq (by integrating the probability 
density P oc e~^^/^ for all values of Hq), and the use of the described independent priors, 
yields the confidence intervals considered below. 

Strong degeneracies exist in the effect of w{z) and on the expansion history, 
requiring independent constraints to make significant progress. Here we consider the use of 
one or more of the following constraints. 

(1) The SDSS luminous red galaxy, baryon acoustic oscillation (BAO) distance parameter to 
z = 0.35: A = QiiE{z)-^l^[\ , where z = 0.35 and A = 0.469(n/0.98)-°-35 ^ q.OI? 
from Eisenstein et al. (2005). The 3-year Wilkinson Microwave Anisotropy Probe (WMAP) 
results yield n = 0.95 (Spergel et al. 2006). 

(2) The present local mass density: VLm = 0.28 ± 0.04, a consensus value when combining 
large-scale structure (LSS) measurements from the 2dF and SDSS of fluh (Tegmark et al. 
2004; Cole et al. 2005) and the value of Hq from EST Cepheids (Preedman et al. 2001; 
Riess et al. 2005). 

(3) The distance to last scattering, z = 1089: In the i/o-independent form (Bond, Efstathiou, 
& Tegmark 1997), Rcmb = Io^^^ ^ ^om the 3-year integrated WMAP analysis 
(Spergel et al. 2006), updated by Wang & Mukherjee (2006) to be 1.70 ± 0.03 independent 
of the dark-energy model. 

(4) The distance ratio from z = 0.35 to z = 1089: Measured by the SDSS BAOs from 
Eisenstein et al. (2005), i?o.35 = ^ps) ^'^'^ [/o ■''(^^'^/^(O-SS))]^^/^)/ /^loss rf^/E(1089) = 
0.0979 ± 0.0036. 

Additional cosmological constraints exist, but in general provide less leverage and may 
be less robust than the above. 

Unfortunately, constraints anchored at 2; = 1089 such as Rcmb and i?o.35 in (3) and (4) 
above require careful consideration when using continuous descriptions of w{z). Specifically, 
we must consider how to evaluate integrals depending on w{z) in a region where we have 
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no data and little intuition, i.e., between the highest- redshift SN la at z 1.8 and the 
surface of last scattering at z = 1089. This span contributes significantly to Rcmb (and 
even more to i?o.35)i accounting for approximately two-thirds of the total distance to last 
scattering. Constraints resulting from a naive leap across this present "cosmological desert" 
are likely to be unjustifiably strong, and extremely sensitive to our assumptions about the 
early behavior of dark energy or on our chosen, parametric form for w{z) (as we will see). 

To evaluate w{z) as well as its sensitivity to assumptions about its high-rcdshift 
behavior, we will make use of three different priors. Our "weak" prior (the most 
conservative) will utilize only VLm and A which arc determined at low rcdshifts, making no 
attempt to guess the behavior of dark energy in the span 1.8 < 2; < 1089. Our "strong" 
constraint integrates the expressions containing w{z) in Rcmb and -R0.35 between z — Q and 
the highest-redshift SN [z ^ 1.8) and beyond; it assumes that the influence of dark energy 
on H{z) is minimal by relaxing the equation-of-state parameter to w — —1 ai z > 1.8. This 
"strong" prior strikes a balance between making maximal use of all information, but only 
broadly guesses at the importance of dark energy where there arc no discrete data (i.e., by 
assuming that the high- redshift Universe is fully dark-matter dominated). Our "strongest" 
prior is the naive extension of w{z) from z = 1.8 to ^ = 1089 strictly along its empirical or 
parametric description. 

The minimum complexity required to detect time variation in dark energy is to add 
a second parameter to measure a change in the equation-of-state parameter with redshift. 
The expansion w{z) = Wq + w'z, where w' = ^\z=o was proposed by Cooray & Huterer 
(1999), OOOfirst used by Di Pietro & Claeskens (2003) and later by R04. However, in this 
form, w{z) diverges, making it unsuitable at high redshift or requiring guesses as to the 
form of its graceful exit (Wang & Mukherjee 2006). Chevaher & Polarski (2001) and Linder 
(2003) suggest w{z) — Wo + Waz/{1 + z), which solves the divergence problem but at a cost 
of demanding stiffer behavior for w{z) than we may presume a priori. As we shall see, 
invoking the {wq^Wo) parameterization can act like an additional, strong prior itself . 

In Figure 9 we show the constraints on wq — Wa plane for the three priors for the Gold 
sample. Here we focus on what is learned from the inclusion of the high-redshift SNe la 
from HST, both the new objects and those with improved cahbration from R04. 

The addition of SNe la at 2; > 1 with the weak prior provides valuable leverage in 
the Wq — Wa space, well in excess of their fractional contribution to the sample. Citing a 
popular metric, the area of the 95% confidence interval (Huterer & Turner 2001, Kolb et al. 

2006; Albrecht & Bernstein 2006), wc note a reduction by 40% or an increase in the figure 
of merit (inverse area) by a factor of 1.7 by the inclusion of the Cold-quality SNe la from 
HST with the weak prior. The reason is readily apparent; because the prior contains no 
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information regarding the evolution of dark energy, what httle is gleaned comes exclusively 
from the discrete data. The highest-redshift SNe la are critical to breaking the degeneracy 
between wq and Wa affecting their lower-redshift brethren. Constraints measured using the 
parameterization, w(z) = Wq + w'z yield similar improvements, 50% in area, when the prior 
is weak. The best fit is consistent with a cosmological constant at Wq = —1 and Wa = 0. 
However, the overall level of empirical knowledge about dark-energy evolution remains very 
modest. We may conclude that if dark energy evolves, the evolution is not very rapid 
(though it remains difficult to predict the natural level of expected evolution; see Caldwell 
& Linder 2005). 

As our prior becomes more ambitious (maintaining the fiatness prior), the confidence 
intervals rapidly shrink. For the strong and strongest prior, the 95% confidence interval 
is 1.5 and 4.4 times smaller than the weak prior, respectively. For the strongest prior, 
even modest values of Wa (and even more so for w') are strongly excluded because of 
their implication that dark energy would gradually and unhaltingly grow in importance 
with redshift to last scattering. For these priors, a cosmological constant (i.e., wq — —1, 
Wa = 0) is separated from the best fit along the direction of the major axis of the error 
ellipse, lying within the boundary of the joint Icr to 2o" confidence level. However, it appears 
that the gains from the stronger priors come at a cost of reliance on the merit of our prior 
at high redshifts and with diminished regard for the data sampling w{z) at lower redshifts. 

We now explore the sensitivity of the constraints on w{z) to the assumption that it 
can be described by a simple, parametric form. To do so we compare the constraints on 
w{z) derived from the previous w{z) = wq + Waz/{1 + z) to a. higher-order polynomial 
expansion in powers of ln(l -|- z) which assumes httle about w{z). We will construct 
w{z) — X^^Wj(ln(l -I- z)Y, using Eq. (3) and a similar statistic to determine the 
likelihood for the terms Wi to i — A (4th order). The seven-dimensional likelihood 
manifold for {wi^Qm, Hq) is computed using a Monte Carlo Markov Chain (MCMC). The 
MCMC procedure involves randomly choosing values for Wi, Qm, and Hq, computing the 
luminosity distance, evaluating the of the fit to supernova data, and then determining 
whether to accept or reject the set of parameters based on whether is improved (using 
a Metropolis-Hastings algorithm). A set of parameters which is accepted to the chain 
forms a new starting point for the next guess and the process is repeated for a sufficient 
number of steps (convergence can be checked via the method of Dunkley et al. 2005). 
Ultimately, the list (aka chain) of the randomly chosen and accepted parameters forms a 



^^For comparison, a constant w model requires w < —0.7 if dark energy does not evolve to last scattering. 
Prom this we can see that even modest evolution, if unwavering, will become incompatible with a high-redshift 
constraint. This result is very similar to that of Wang & Mukerjee 2006. 
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good approximation of the likelihood distribution (Knox, Christensen, & Skordis 2001). 

The best solution and its uncertainty for w{z) is shown in Figure 10 for the quartic 
polynomial and low-order descriptions oi w{z). As shown, the low-order fit implies much 
greater precision concerning dark energy (with a 95% confidence region which is 7.2 times 
smaller) . 

The higher-order fit suffers in comparison, tolerating non-monotonic and even 
oscillating solutions for w{z) as well as recent changes. The only (natural) limitation to the 
high-order fit comes from direct confrontation with data. Inclusion of the iJST-discovered 
SNe la alone reduces the confidence intervals by a factor of 1.6. Yet, despite the extra 
freedom, the recent equation- of- state parameter remains well constrained. In Figure 11 we 
show the uncertainty as a function of redshift for w{z) for the low-order and high-order 
fits with and without the high-redshift HST data. The difference in implied precision on 
w{z) from the two parameterizations is very similar to the difference seen using a weak or 
strong prior at high redshift. Both a simple parameterization and a strong, high-redshift 
prior greatly restrict the allowed wandering of w{z), though neither can be well-justified. 
In this way (and others discussed in the next section) we can see that a simple dark-energy 
parameterization is equivalent to a strong and unjustified prior on the nature of dark energy. 
The conclusions we draw either from the polynomial form with a strong prior or from the 
simple form with a weak prior are powerfully shaped by the data and are unlikely to be 
completely mislead by a bad assumption. 

An alternative approach to parameterizing w{z) and constraining the parameters is to 
extract discrete, uncorrelated estimates of w{z) as a function of redshift, analogous to the 
uncorrelated estimates of H{z) derived in the last section. Following a method established 
by Huterer & Cooray (2005), we can extract the uncorrelated and model- independent band 
power estimates of w{z) similar to a principal component analysis. 

The evolution of w may be usefully resolved across several redshift bins or steps. Let 
Wi be a constant value of w{z) in bin i. For a given experiment, the Wi will generally be 
correlated with each other because measurements of di constrain redshift integrals of w{z); 
see Eq. (8). We can use a MCMC to obtain the likelihood surface for the variables in the 
vector w = [wi, W2, ■ ■ ■ Wn]- 

Assume for simplicity that the mean value of w from the likelihood surface has already 
been subtracted so that now (wj) = 0. In general, the variables will be correlated, with a 
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correlation matrix 

/ < 

C = (ww^) = ; 

What we're interested in is finding a linear transformation from w into some new vector 
VV such that the correlation matrix for the new vector is diagonal (in other words, we want 
(yViVVj) = for i 7^ j). Since there are an infinite number of matrices that diagonalize C 
there are an arbitrary number of ways to do this. 

The simplest is to decompose C into eigenvalues and eigenvectors, and use the 
eigenvectors to multiply w. If 

C = V^AV, 

where A is diagonal, then in turn 

VCV^ = A. 

— * 

Therefore if we choose the vector W to be equal to Vw, then 

(WWT) = ((Vw) (Vw)^) 

= V(ww^)V^ = VCV^ = A. 

— * 

So in this case W has a diagonal correlation matrix. However, transforming from w 

— * 

to W in this way will often involve both adding and subtracting individual entries of w to 

— * — * 

get a particular entry of W so physically interpreting what W means can be difficult. It's 
actually possible to use a different transformation to get a W that is in some sense "closer" 
to w. A slightly more complicated transformation is to define a matrix T such that 

W = Tw = V'^A-^Vw, 

where A~ 2 is just the diagonal matrix of the reciprocal of the square root of each eigenvalue. 

— * 

With this transformation the covariance matrix for W is 



(WW^) = ((V'^A-Ww) (v^A-iVw) ) 

= V'^A-iv(ww^)V^A-W 
= V'^A-^VCV^A-iv 
= V^A-^AA-^V = I. 

— * 

So with this transformation the correlation matrix for W is just the identity matrix. 
Further, with a matrix of the form V'^A^a V, the entries of T are much more likely to be 
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positive (the eigenvalues are positive, so any negative entries in the eigenvector matrix V 
will tend to multiply together and become positive). 

This is the transformation we used, with one further modification. Explicitly writing 

— * 

out the transformation for a single entry in W 

j 

For a given i, Tij can be thought of as weights for each Wj in the transformation from w 
to W) . We are free to rescale each W) without changing the diagonality of the correlation 
matrix, so we then multiply both sides of the equation above by an amount such that the 
sum of the weights J2j Tij is equal to one. This allows for easy interpretation of the weights 
as a kind of discretized window function. This is, in fact, what we plot as the window 
function for the decorrelated parameters in the paper. 

The resultant values of Wj are thus uncorrelated (in that there is no covariance between 
the errors of any pair) and free to vary independently, subject only to the constraints of 
the data. Each uncorrelated measure of Wj can be tested against the cosmological constant 
expectation of Wj = — 1, independently as well as collectively. Again, we use the weak, 
strong, and strongest priors, the last two including Rcmb and -R0.35 with w{z > 1.8) 
equal to —1 or the value in the highcst-redshift bin, respectively. Wc also used the same 
three rcdshift bins defined by nAz = 40, resulting in three independent measures of w{z) 
anchored at approximately z=0.25, 0.70, and 1.35 (with boundaries of z = 0, 0.45, 0.935, and 
1.8) for the terms W0.25, Vyo.7o, and W1.35, respectively. Thus there are five free parameters 
in the MCMC chain corresponding to the 3 values for Wj, Qm, and Hq. The likehhoods 
and confidence intervals for these terms are shown in Figure 12, 13, and 14 (respectively) 
and given in Table 5. The rows of T are represented as the window functions in Figures 12, 

13, and 14. As expected, the lowest-redshift measure of w{z), Wo. 25, is derived primarily 
from the lowest-redshift bin which provides over 90% of the weight. The highest-redshift 
measure, VV1.35, is mostly derived from the highest-redshift bin which contributes two-thirds 
of the weight of the measurement, but also relies on the lower-redshift bins to decouple 
high-redshift measurements of di from the low-redshift behavior of w{z). 

For the two lowest-redshift bins, the likelihood distributions are close to Gaussian, but 
for W1.35 the distribution can be quite skewed, requiring an explanation of the definition of 
our confidence intervals and reported values of . In Table 5 and Figures 12, 13, and 

14, we define a +la and —la region to be the boundaries at which the likehhood falls to 
0.6065 of the peak on either side (i.e., as for a Gaussian), a definition which neglects the 
non-Gaussian tails to provide the frequently sought error bar. Likewise, we define a 2a 
region as the boundaries of equal hkelihood which contain 95% of the hkelihood. We also 
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report the values of the peak of the hkehhood for in the table. 

For the weak prior, the three measures (>Vo.25, VV0.70, and W1.35) are all consistent 
with w = — 1 at or near the 68% confidence interval. The lowest- redshift bin remains well 
constrained with W0.25 = —1.06 ± 0.10 despite the additional freedom at higher redshifts. 
The next bin is modestly well constrained at W0.70 = —0.46 ± 0.46, though the hkelihood in 
the highest-redshift bin has only begun to appear at all localized. Although the maximum 
likelihood values of w{z) are monotonic in redshift, any trend is not significant. The 
constraints are weakest at high redshift where the data are sparse and discrete values of 
w{z) become harder to isolate from integrated constraints. Figure 15 shows the impact of 
adding SNe la to the sample at the highest redshifts from HST. The first published sampled 
from R04 began the process of locahzing the highest redshift bin and markedly improved 
the intermediate bin, with no impact on the nearest bin. For the weakest prior which lacks 
any complementary constraints at high redshift, the increased sample has provided only 
modest gains. 

For the strong prior, the nearby bin is unchanged, the uncertainty in the intermediate 
bin tightens by 25%, and the highest-redshift bin becomes significantly more peaked with a 
factor of two reduction in uncertainty. A modest tension exists with w — —1 in the middle 
bin (less than 2a) (with no apparent trend in redshift). However, most of the data in this 
intermediate bin come from previous and on-going ground-based surveys and is rapidly 
growing in size. Thus we can expect this bin to soon come into better focus. Interestingly, 
the full sample considered here will soon be significantly augmented by the addition of 60 
new SNe la with high-quality data from the first ESSENCE dataset of Wood-Vasey et 
al. (2007) which taken together with the past data, we note, results in a reduction of the 
difference with w — —1 at these redshifts to ~ 1 a. 

For the strong prior, the addition of the high-redshift HST data shown in Figure 15 
provides strong gains in the precision of W1.35, especially with the sample published here 
for which we find W1.35 = — O.Sl? q, marking significant inroads into the measurement of the 
early-time behavior of dark energy. Indeed, this measurement represents the most distinctive 
contribution of the iJST-discovered SNe la. Combined with the low-redshift measurement, 
Wo.25, the change Wo. 25 — VV1.35 = — 0.3l5;o over one unit of redshift, disfavoring rapid 
changes, i.e., \dw/dz\ » 1. Indeed, without the HST data, any measurement of the 
same Wi,35 is not meaningful because the sample would contain no SNe at z ~ 1.35. 
The gain from the weak to the strong prior results from the complementary nature of 
the contributions of the SN data, the BAO measurement, and the WMAP measurement 
moderated by the high-redshift integration. 

For the strongest prior, the tension within the highest-redshift bin has grown between 
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the SNe contained within and the extrapolation of the imphed behavior of dark energy 
continuously to z = 1089. Any growing mode or non-accclcrating type behavior for dark 
energy (w > —1/3) is curtailed. While this plot looks convincing, it is worth noting that 
non-accelerating modes were mildly favored by the weak prior. The embedded assumption 
of the way dark energy changes with redshift dominates over the information added by 
considering the data on w{z). Because this prior relies on the correctness of guesses, we 
conclude that the apparent gain in information is at best risky and at worst, illusionary. 
Thus, we discourage the use of other parameterizations of dynamic dark energy in which 
the behavior of dark energy between z — 1.8 and z — 1089 remains important and a 
determining factor. 

Not surprisingly, the addition of the SNLS sample to the full sample provides the most 
improvement in the lowest-redshift bin, about 10% in a^, but little gain in the intermediate 
bin (3%) and none in the highest redshift bin. 

The likelihood distribution for Vyi.35 also empirically characterizes the most basic 
property of dark energy at z > 1, the way it gravitates. Using the strong prior, the 
hkelihood that high-redshift dark energy has negative pressure, i.e., W1.35 < 0, is 97.6% (or 
99% for the MLCS2k2 fits to the Astier et al. light curves). (We note that the likehhood 
contained beyond a single boundary such a w < cannot be directly inferred from our 
non-gaussian confidence regions because they have been defined using two boundaries and 
the requirement that the boundary likelihoods match.) Thus, the defining property of 
dark energy appears to be intact even during the epoch of matter domination. We can 
say with greater confidence (99.99%) that if a decaying scalar field is responsible for this 
high-redshift dark energy, its energy is not primarily kinetic {w = 1). With less confidence, 
93.4%, we can say that early dark energy provided repulsive gravity {w < —1/3). The 
nature of dark energy at high redshift {z > 1) may be of particular interest because 
classes of tracker models predict rapid changes in w{z) near epochs when one energy 
density (e.g., radiation, dark matter, or dark energy) becomes dominant over another. For 
example, Albrecht and Skordis (2000) describe a form of potential which gravitates as the 
dominant mass-energy component at a given epoch. Simply put, this model would predict 
that the equation of state of the dark energy would mimic matter {w = 0) during the 
decelerating, matter-dominated phase (with a brief transition "ripple" during which w > 0) 
and the aforementioned likelihood for W1.35 > would appear inconsistent with this model. 



^°The strongest prior would appear to provide even greater confidence for this statement. However, as we 
discussed, this prior seems unjustifiably strong and its use in this context is somewhat circular (i.e., guessing, 
and then inferring, the nature of dark energy at high redshift). In contrast, the strong prior brings the most 
independent information to bear on the early-time behavior of dark energy. 
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However, future observational and theoretical work should better define if variants of this 
idea remain viable. 

Finally, we may consider whether three additional parameters to describe w{z) are 
actually needed to improve upon a flat, A-cold-dark-matter (ACDM) model fit to the data 
(i.e., the "concordance" value of — 0.29, JIa = 0.71). The residuals from this fit are 
shown in Figure 17 and have a dispersion for all data of 0.21 mag. To determine the need 
for complex forms of 'u;(^) we can calculate the improvement to the fit, 

xlff = -A(21n£) = 21n£(w = -1) - 2lnC{wi = Wi), (9) 

with i additional free parameters. For the weak, strong, and strongest priors we find an 
improvement of xlff = 4, 5.5, and 5.5, respectively, for the three additional degrees of 
freedom, in no case requiring the additional complexity in dark energy (improvements 
of > 14 would be noteworthy). Likewise, there is no improvement at all for the Akaike 
Information Criterion (i.e., AAIC=Ax^ — 2i; Liddle et al. 2004) with changes of —2, —0.5, 
and —0.5, respectively which fail to overcome the penalty of increased complexity in the 
model. 



4. Discussion 

We have presented a sample of SNe la at 2; > 1.0, found and followed with HST, more 
than double the size of those in R04. We have also performed a set of limited cosmological 
analyses focusing on tests and models we expect to be most sensitive to the additional data. 
Given the set of independent experimental constraints and the panoply of dark-energy 
models, we expect that the most interesting tests and constraints will be discovered by 
future work. We now discuss caveats and tests of this dataset. 



4.1. Systematics 

Past research has sought to quantify and characterize possible systematic errors in 
supernova-based cosmological inference (e.g., Riess et al. 1998; Perlmutter et al. 1999; 
Knop et al. 2003; Tonry et al. 2003; Riess et al. 2004; Filippenko 2005; Astier et al. 2006). 
Here we consider systematics which are more sensitive to the uniquely high redshifts of our 
sample or otherwise merit new consideration. 

Strong or weak lensing may skew the distribution of magnitudes but will not bias the 
mean (Holz 1998). Extensive analyses of the expected lensing along the line-of-sight for the 
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SNe found in the GOODS fields in R04 by Jonsson et al. (2006) show that the distribution 
of magnifications matches an unbiased sample having a mean consistent with unity and 
that attempting to correct for this effect has negligible impact on the derived cosmological 
parameters. This is not surprising, as the dispersion of lensing-induced magnitudes is less 
than 0.05 mag, far less than the intrinsic dispersion. The same was concluded in R04 using 
the methods of Benitez et al. (2002), though with less accuracy than in Jonsson et al. 
(2006). 

As high-rcdshift SN la surveys continue, sample sizes increase, and data collection 
methods improve, it may become sensible to discard older samples in an effort to reduce 
systematic errors in the cosmological sample associated with inferior calibration (e.g., data 
from Riess, et al. 1998, Perlmutter et al. 1999, Tonry et al. 2003, and Knop et al. 2003). 
At present we do not believe this is warranted, especially since it is not yet possible to 
similarly replace the low-redshift sample whose data was collected with similar techniques 
as the older high-redshift data. Nevertheless, it may be instructive to determine the impact 
on the determination of the nature of dark energy using at intermediate redshifts only the 
more modern and homogeneous data from Astier et al. (2006) and with the same HST data 
at high redshift and low-redshift sample found here. From this more limited sample we find 
no significant change for the lowest redshift measure of the equation of state, W0.25, and a 
30% decrease in its precision. For the intermediate redshift bin, Wo.to, we find values more 
consistent with a cosmological constant with — l.ll+oJa and — 0.76jlQ;|g , for the weak and 
strong prior, respectively, representing a 50% to 100% decrease in the precision. For the 
highest redshift measure, W1.35, we find values less consistent from a cosmological constant, 
with 3.391^9! and -0.23l[j:?^ , for the weak and strong prior, respectively, representing little 
change in the mean precision. The two most discrepant measurements from w{z) = —1 are 
with the full datasct and the strong prior for VVo.7o=-0.15 (-lo'=-0.41, -2 a =-0.89) and for 
the limited datasct and the weak prior for which Wi.35=3.39 (-1 o"=-2.91, -2 a=-23.1). Both 
measurements are inconsistent at more than the 1 sigma level and less than the 2 sigma level 
with a cosmological constant, an insignificant difference for 3 independent measurements 
of w{z). We also note that an increased sample of 60 new SNe la at 0.3 < 2; < 0.7 
from ESSENCE (Wood-Vasey et al. 2007) further reduces the apparent difference from a 
cosmological constant from the full-sample for the intermediate redshift bin. 

Evolution of SN explosions remains the most pernicious systematic source of uncertainty 
and is challenging to quantify (see, however, Riess & Livio 2006 for a suggestion for future 
observations). The two most direct lines of evidence that evolution with redshift, is limited 
to less than ~10% in luminosity comes from studies of SN la distances across host-galaxy 
(characteristic) ages (Sullivan et al. 2003; Riess et al. 1998; Jha 2002; Jha, Riess, & 
Kirshner 2006, Gallagher et al. 2005) and the consistency of SN la distances with the 
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expected deceleration at z > 1 (R04). 

An alternative to evolution proportional to redshift would be evolution that is 
proportional to look-back time (Wright 2002). However, this seems unhkely, as the 
evolution of cosmic properties to which SNe la may respond is more closely correlated with 
redshift than with time. Two natural examples are the evolution of cosmic metallicity as 
seen in damped Lya systems (Kulkarni ct al. 2005) and the minimum stellar mass of a star 
which can turn off the main sequence and donate the SN la white dwarf. Combined with 
the necessary delay required to produce SNe la (> 1 Gyr), the gradient in these physical 
parameters lies primarily with redshift, not time (Riess & Livio 2006). 

Recent work has questioned the validity of the mean Galactic reddening-law parameter, 
Ry = Ay / Eb^v = 3.1, for use in determining the extinction in extragalactic SNe. 
Estimates of the typical value of Ry in distant galaxies hosting SNe la tend to prefer values 
somewhat less than 3.1 (Riess, Press, & Kirshner 1995; Jha 2002; Jha, Riess, & Kirshner 
2006; Philhps et al. 1999), with values as low as Ry = 2.3 by Wang et al. (2006). An 
apparent correlation exists between Ry and Ay: SNe with greater extinction prefer lower 
values of Ry, and both of the well-studied SNe la with Ay > 2 mag (SN 1999cl and SN 
1996ai) yield Ry < 2.0 (Jha, Riess, & Kirshner 2006). Wang et al. (2006) have suggested, 
as a natural explanation for reduced values of Ry, the occasional dominance of scattering 
over absorption caused by dust clouds in the circumstellar environment. Evidence does 
exist for the occasional and significant presence of circumstellar dust from at least 3 SNe la 
exhibiting light echos (SN 1995E, Quinn et al. 2006; SN 1998bu, Gappcllaro et al. 2001; SN 
1991T, Schmidt et al. 1994) and a few possible SNe la with circumstellar hydrogen seen at 
late times (e.g., SN 2002ic; Hamuy et al. 2003). In the presence of local dust, Wang et al. 
(2006) show that 1.8 < Ry < 3.1 for normal. Milky- Way-type dust. 

However, past determinations of the mean Ry for distant SNe la are strongly biased 
toward the few objects with large extinction and lower Ry, which have greater weight in 
these analyses. For example, Jha, Riess, & Kirshner (2006) analyzed the 33 nearby SNe la 
with the highest extinction {Ay > 0.5 mag) and found a median Ry — 2.8, well in excess 
of the weighted mean value of 1.8 (with half of this difference explained by the two SNe la 
with Ay > 2 mag). Because our cosmological analysis excludes (from the Gold sample) all 
SNe with Ay > 0.5 mag, we would expect the characteristic Ry in our sample to be better 
represented by the median, or perhaps closer to Galactic for low-extinction events. Further 
analyses of large samples of low-extinction events may yield greater insight on this issue. 

Analyses which seek to correlate SN color with luminosity, such as those by Tripp & 
Branch (1999), Guy et al. (2005), and Astier et al. (2006), combine both the reddening-law 
parameter Ry and intrinsic luminosity-color relation parameter into a single relation. 
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If these independent, physical relations have the same dependence between color and 
brightness, or if the relative contributions of dust and populations to the relation do not 
change with redshift, then combining these effects into one should succeed. However, if 
one of these assumptions fails significantly (e.g., there is more or less dust on average at 
different redshifts), important systematic errors in the cosmological inference may result. 
In any case, the effective value of Ry from these analyses should be an average of the value 
from dust and the physics of SN explosions. The MLCS2k2 templates suggest that the 
intrinsic relation from normal SNe la would be M/{B — V)!v 1.8, reducing the effective 
Rv below the typical value for the Galaxy as found by Tripp & Branch (1999), Guy et al. 
(2005), and Astier et al. (2006). 

To determine the sensitivity of our cosmological analyses to a reduced value of Ry, we 
refit all of the distances using Ry — 2.5 and 2.0 and compared the results to our preceding 
distances determined with Ry — 3.1. We found our mean distance to increase by 0.026 
mag and 0.042 mag for Ry = 2.5 and Ry = 2.0, respectively. However, we found the 
mean variation of this change to be negligible across redshift with less than a 0.01 mag 
difference between the bins used in Figure 7. The distance determinations to objects with 
large extinction {Ay > 1 mag) depend more critically on Ry, with an increase of 0.10 and 
0.35 mag for SN 2002kc {Ay = 1.35 mag for Ry = 3.1) for Ry=2.5 and 2.0, respectively 
Such sensitivity is exactly the reason it is unwise to make use of well-reddened objects in 
cosmological analyses and why we limited our Gold sample to Ay < 0.5 mag. 

Wood-Vascy et al (2007) has shown that a selection bias exists in the distribution 
of observed extinctions and luminosities of SNe la found in the ESSENCE Survey for 
supernovae discovered near the redshift hmit of the survey. This bias can be corrected using 
Monte Carlo simulations of the survey selection. The bias becomes significant for SNe la 
which can only be found at a signal-to- noise ratio of less than 15 due to their apparent 
faintness and the survey's hmitations. For our HST SN survey, this hmit occurs for SNe la 
at ^ > 1.4 and thus such selection bias has little affect on our sample, the great majority of 
which is at 2; < 1.4. The reason for the absence of SNe la at the redshift limit of our survey 
{z ~ 1.7) has been attributed to an apparent 2 to 3 Gyr delay between the formation of 
stars and the production of SNe la (Strolger et al. 2004). 

At present, none of the known, well-studied sources of systematic error rivals the 
statistical errors presented here. 

Now we consider systematic errors whose definition includes our present inability 
to discover them and to correct for them. Here we will quantify the translation of such 
systematic errors to our data analysis and their propagated effects on the cosmology. A 
simple form for general systematic errors in SN distances to take is a correlation with 
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redshift. Such a model provides an adequate (if naive) description of zeropoint errors tied 
to fixed passbands, errors in K-corrections, random errors from evolution (tied to look-back 
time), etc. This model is good at addressing zeropoint errors between small ground-based 
telescopes used to collect the z < 0.1 sample, errors from large ground-based programs 
which collected the sample at 0.1 < z < 0.8, and from HST at 2; > 1. It is similar to the 
model used by the Supernova Acceleration Probe (SNAP) collaboration (Aldering 2005). 

We define the covariance between two SNe to be < o-^,iO-^,2 >= A^cxp[— |2;i — Z2\/zt], 
where A provides the correlation (in mag) for two SNe at the same rcdshift and zt provides 
the decay length of the correlation in redshift. For the expected precision of calibration of 
HST and past datasets, we examine a model with {A, Zt) = (0.05,0.05), and (0.05,0.10) 
against the error-free case (0.0,0.0). We consider (0.05,0.10) to be a rough estimate of the 
covariance of the SN dataset, in which SNe at hke redshifts suffer the same 0.05 mag errors 
related to cahbration, those transformed to adjacent rest-frame bandpasses {zt > 0.20) 
share 0.02 mag of systematic error, and the entire Hubble flow sample {zt = 0.1) shares an 
0.03 mag error in common (and is not improved beyond the inclusion of ~ 100 SNe la). 
For the more optimistic model, (0.05,0.05), the correlation length of these errors is halved 
in redshift. 

For the Gold sample of SNe la presented here, we use the redshifts and the distance 
errors (but not the distance measurements) to compute the Fisher Information Matrix 
combined with a typical prior on ctq^^ = 0.03. To study these systematic errors we include 
a covariance matrix with off-diagonal terms given by our model and parameters. We then 
compute the expected 68% confidence level interval for the parameter space Wq — w' for 
the 3 cases. We find that our best guess at covariance, (0.05,0.10) broadens the interval 
by 27% whereas the more optimistic (0.05,0.05) broadens the interval by 17%. Thus we 
conclude that our present results remain dominated by statistical errors. However, as the 
sample grows, for this balance to remain true, ongoing improvements in calibration must be 
realized. 



5. Summary and Conclusions 

(1) We present 21 new if.ST-discovered SNe la and an improved calibration of the 
previous sample from R04. Together this sample contains 23 SNe la at 2; > 1, extending the 
Hubble diagram over 10 Gyr. 

(2) We derive uncorrelated, model-independent estimates of H{z) which well-dehneate 
current acceleration and preceding deceleration. The iJST-discovered SNe la measure 
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H{z > 1) to slightly better than 20% precision. 

(3) The full if-ST-discovered SN la sample, presented here, provides a factor of 
two improvement over our present ability to constrain simple parameterizations of the 
equation-of-state parameter of dark energy (w) and its evolution. 

(4) Stronger priors and tighter constraints on the preferred cosmological model can be 
extracted from independent measurements tied to the surface of last scattering, but the 
use of these requires assumptions about the behavior of dark energy across a wide range 
of redshift (1.8 < z < 1089). The strongest of these priors, hke the simplest dark energy 
parameterizations, appears unjustified in the presence of our current ignorance about dark 
energy. Assuming the effect of dark energy at 2; > 1.8 is minimal, we derive meaningful 
constraints on the early properties of dark energy: 'w{z > 1) = — 0.8l°;o and 'w{z > 1) < 0, 
i.e., negative pressure, at 98% confidence. 

(5) At present, we find that the use of additional parameters to describe w{z) does not 
provide a statistically significant improvement to the fit of the redshift-magnitude relation 
over the use of a simple cosmological constant. 

(6) An analysis of the z > 1 sample-averaged spectrum shows it to be consistent with 
the mean spectrum of SNe la over the last 10 Gyr, failing to reveal direct evidence for SN la 
evolution. 
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AURA, Inc., under NASA contract NAS 5-26555. Some of the data presented herein were 
obtained with the W. M. Keck Observatory, which is operated as a scientific partnership 
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Table 4. Comparison of Gold Set Data to Models, cz > 7000 km s~^ 



Model x^(for 184 SNe la) 



= 0.29, = 0.71 150« 

JIm = 1-00, Oa = 0.00 285« 

= 0.00, Qa = 0.00 164« 

High-redshift gray dust (with = 1-00, Qa = 0.00) 344^ 

Replenishing dust (with Qm = 1-00, Qa = 0.00) 150^ 

Dimming oc z (with Vtu = 1-00, Qa = 0.00) 266'' 



^best after marginalizing over Hq 

^hesi for best Hq 



-51 - 



Table 5. Likelihood Regions For Wz 
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la 2 a 








Prior=Weak, Sample=All Gold 
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Wo.70 
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Wl.35 


0.59 


-2.62 to 3.03 -16.6 to 6.15 






Prior=Weak, Sample=Gold minus HST 


Wo.25 
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0.11 
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10.77 


1.86 to 18.55 -20.1 to 27.92 






Prior=Strong, Sample=All Gold 
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-0.15 
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Prior=Strongest, Sample=Gold minus HST 
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-1.03 


-1.13 to -0.94 -1.24 to -0.85 
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0.24 


-0.17 to 0.64 -0.70 to 1.06 
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-1.89 


-5.50 to -0.80 -18.0 to -0.34 




Prior=Strong, Samplc=All Gold with MLCS2k2 Fits to SNLS SNe 


Wo.25 


-1.05 


-1.14 to -0.94 -1.26 to -0.84 


Wo.70 


-0.09 


-0.45 to 0.23 -0.91 to 0.56 


Wi.35 


-1.01 


-2.23 to -0.26 -15.8 to 0.37 
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Table 6. HST-discovered Sample*; Distance Scale of Riess et al. 2004 
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Figure Captions 

Figure 1: Discovery- image sections from ACS FS50LP images around each SN. Panels 
on the left and middle show the discovery epoch and the preceding (template) epoch, 
respectively. The panels on the right show the results of the subtraction (discovery epoch 
minus template). Arrows indicate position of the SNe. Image scales and orientations are 
given. 

Figure 2: Color images of hosts with site of SN indicated. Filters F850LP, F6061^, 
and FA35W correspond to red, green, and blue, respectively. 

Figure 3: Corrections to single-model PSF fitting photometry in the F850LP band for 
SNe la due to the red halo. At red wavelengths, the PSF model and aperture corrections for 
ACS WFC strongly depend on the monochromatic wavelength of incident hght as measured 
from narrow-band filters. Thus the use of a single PSF model or aperture correction 
to measure photometry in a wide passband necessitates the tisc of a correction derived 
synthetically by weighting the monochromatic photometric variations by the appropriate 
SN SED. In the upper panel we show this correction for F850LP as a function of SN la age 
for different redshifts: 0.75 (asterisks), 1.00 (diamonds), 1.30 (squares), 1.50 (Xs), and 1.70 
(triangles). In the lower panel we show the correction as a function of redshift for a SN la 
at maximum light. 

Figure 4: Multi-color light curves of SNe la. For each SN la, multi-color photometry 
transferred to rest-frame passbands is plotted. The individual, best-fit MLCS2k2 model is 
shown as a solid line, with a ±la model uncertainty, derived from the model covariance 
matrix, above and below the best fit. 

Figure 5: Identification spectra (in fx) of 12 of the new i/ST-discovered high-rcdshift 
SNe la all shown in the rest frame. Classification features are analyzed in §3. The data are 
compared to nearby SN la spectra of a similar age as determined by the light curves (see 
Table 3). 

Figure 6: MLCS2k2 SN la Hubble diagram. SNe la from ground-based discoveries 
in the Gold sample are shown as diamonds, if^T-discovered SNe la are shown as filled 
symbols. Overplotted is the best fit for a flat cosmology: VIm = 0.27, Ha ~ 0.73. Inset: 
Residual Hubble diagram and models after subtracting empty Universe model. The Gold 
sample is binned in equal spans of nAz — 6 where n is the number of SNe in a bin and Az 
is the redshift range of the bin. 

Figure 7: Uncorrelated estimates of the expansion history. Following the method of 
Wang & Tegmark (2005) we derive 3, 4, or 5 independent measurements of H{z) from the 
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Gold sample using n/S.z — 40, 20, and 15, respectively. The bottom panel shows the derived 
quantity a versus redshift. In this plane a positive or negative sign of the slope of the data 
indicates deceleration or acceleration of the expansion, respectively 

Figure 8: Same as upper panel of Figure 7 comparing the improvement to the 
highest-redshift measure of H{z) due only to the newest HST data, i.e. since R04. 

Figure 9: Joint confidence intervals derived from SN samples for a two-parameter 
model of the equation-of-state parameter of dark energy, w{z) — Wa + Waz/{1 + z). For 
each panel, constraints from a SN sample are combined with the indicated prior to yield 
the indicated confidence intervals. The position of a cosmological constant (—1, 0) is shown 
as a filled symbol. 

Figure 10: Best solution and uncertainty for a quartic polynomial fit (blue), and simple 
— Wa parameterization using a strong prior (red). As seen, the simple parameterization 
highly constrains the behavior of w{z) as compared to the polynomial. The greater 
constraint on w{z) implied by the Wq — Wa parameterization derives from implicit and 
unjustified priors on dark energy: that its evolution is monotonic, linear, and most 
important at low redshifts. 

Figure 11: As in Figure 10, uncertainty for w{z) for wq — Wa and quartic polynomial 
parameterizations with and without high- redshift HST data. 

Figure 12: Measurement of 3 uncorrelated components of w{z) using the Gold sample 
of SNe la and the weak prior. Following the method of Huterer and Cooray (2005) we 
derived measurements of w{z) in the same 3 redshift bins used in Figure 7 {nAz=4:0). 
Using their covariance matrix we derived new, uncorrelated components of w{z) shown in 
the upper panel with window functions given in the lower right panel and likelihoods given 
in the lower left panel. 

Figure 13: Same as in Figure 12 but using the strong prior. 

Figure 14: Same as in Figure 12 but using the strongest prior. 

Figure 15: Difference between Figure 12,13 and 14 with and without the HST-discovered 
SNe la for the weak, strong and strongest prior. However, without the HST data, any 
measurement of the highest-redshift value, W1.35, is not very meaningful because the sample 
would contain no SNe at 2; ~ 1.35. 

Figure 16: Average spectrum derived from HST ACS grism spectra of thirteen SNe 
la z > 1. The high redshift average and dispersion (mean z — 1.1) shown in thick gray 
compares well to the low redshift average (over the 10 days following maximum) and 



-78- 



dispersion (day-to-day, for the 10 days post maximum) shown as the heavy hne with error 
bars. 

Figure 17: Distance difference in magnitudes for all Gold SNe {cz > 7000 km 
s"^) between the measured distance and that predicted for the concordance cosmology 
{Qm = 0.29, = 0.71) with w — —1. As discussed in the text, increased complexity in the 
description of w{z) is not presently justified by the improvement in the fit. 
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Table 6 — Continued 



SN 


z 


MS 


a*- 


host Ay 


A 


sample 





*The full sample used here for the cosmological analyses consists of the union of Table 5 from Riess et al. 
2004, this table (with the new distances to SNe from 2002 to 2003 replacing those in Riess et al. 2004), and 
the SNLS sample from Astier et al. (2006) listed in §2.3 of this paper. The full sample is also available at 
http://braeburn.pha.jhu. edu/^ariess/ROe or upon request to ariess@stsci.edu. 

^Redshift and velocity error and intrinsic SN la dispersion of 0.08 mag already included. 

^Distance normalization is arbitrary; see Appendix. 



